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then  used  to  exanine  the  theoretical  basis  for  the  disoriainants  between 
earthquakes  and  underground  nuclear  e^losions  based  on  the  relative 
aaplitudes  of  P  end  S  WKvea.  The  ratio  of  Sn  to  Pn  aaplitude  looks 
proaising  as  a  high  frequency  discrininant.  However,  the  ratio  of  Lg  to 
Pn  ap^litudes  is  not  as  useful  because  the  ocaplex  nature  of  the 
propagation  characteristics  of  Igt.  r”  j 

An  ispleaentation  of  the  operator  aethod  by  aeans  of  coupled  nodes  has 
been  used  to  exanine  the  effect  of  heterogeneous  crustal  structures  on  Lg. 
In  this  technique  the  local  seisaio  wavefield  in  the  real  nediun  is 
expresaed  as  a  oopbinatian  of  the  pedal  eigenfunotions  of  a  stratified 
reference  structure.  Departures  of  the  seismic  pre^rties  in  the  real 
nediun  from  those  in  the  reference  lead  to  coupling  between  the  amplitude 
coefficients  in  the  n^al  exi^sion.  The  evolution  of  these  nodal 
weighting  functions  with  horizontal  position  are  described  by  a  coupled 
set  of  ordinary  differential  equations.  This  approach  provides  a 
oaloulation  soheaie  for  studying  guided  wave  propagation  over  extended 
distances,  at  frequencies  of  1  Hz  and  above.  The  heterogeneity  models 
which  have  been  used  are  two-dimensicxial  and  calculations  are  carried  out 
for  one  frequonoy  at  a  time. 


One  of  the  major  effects  of  crustal  heterogeneity  is  to  introduce  the 
possibility  of  changing  the  character  of  the  main  amplitude  peak  in  the  Lg 
wavetrain  by  shifting  energy  between  different  group  velocity  components. 
As  a  result,  an  effeotive  measure  of  the  oraergy  content  of  the  Lg  waves 
will  be  to  consider  the  integrated  energy  al(xig  the  trace  between  group 
velocities  of  3.6  and  3.3  km/'s.  The  effects  of  heterogeneity  vary  between 
different  parts  of  the  Lg  wave  train  and  can  be  visualised  by  using 
synthetic  seismograms  computed  with  a  narrow  band  response  in  frequency, 
with  allowance  for  interaction  between  the  modes.  The  representation  of 
the  wavefield  in  terms  of  modal  contributions  allows  a  detailed  analysis 
in  terms  of  the  group  velocity  components. 
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RESBARCH  GBJBCTIVBS 

The  ain  of  this  work  is  to  develop  techniques  which  can  be  used  to 
describe  the  propagation  of  regicxial  seisnic  phases  in  three-dijnensionally 
heterogeneous  media,  in  order  to  improve  the  understanding  of  the  nature 
of  the  phases  and  the  way  in  which' the  characteristics  of  the  seisnic 
source  can  be  modified  by  propagation  to  the  receiver.  Such  information 
will  be  valuable  in  assessing  the  effects  of  geological  structure  on  the 
behaviour  of  potential  discriminants  between  earthquakes  and  underground 
nuclear  explosions. 

Host  of  the  regional  phases  of  interest  for  nuclear  discrimination 
problems  are  observed  after  propagation  through  considerable  distances. 
These  phases  travel  through  the  crust  and  uppermost  mantle,  which  from  a 
wide  variety  of  studies  are  known  to  be  regions  of  ctmsiderable  horizontal 
variability  of  properties. 

We  have  therefore  set  out  to  find  theoretical  descriptions  of  the 
propagation  of  the  regional  phases  which  enable  us  to  find  general  results 
on  the  relative  amplitudes  of  different  phases.  In  addition  we  have 
undertaken  a  computational  study  of  the  interaction  of  the  Lg  phase  with 
specific  models  of  heterogeneity 


RBSBARCH  STATUS  AND  AdUKVUUflri'S 

In  general,  observations  of  regional  seismic  phases  are  made  at  some 
hundreds  of  kilometres  from  the  scarce,  so  that  the  waves  can  have  made 
quite  complex  interactions  with  three-dimensional  heterogeneity  in  the 
crust  and  uppermost  mantle.  The  descriptions  of  the  propagation  process  is 
a  difficult  mathaaatical  task  but  can  be  simplified  by  introducing  the 
concept  of  propagation  operators  (Rennett  19^).  This  approach  allows  the 
sequence  of  physical  processes  from  the  generation  of  the  seismic  waves  at 
the  source  to  reception  at  the  sensors  to  be  represented  as  the  action  of 
a  set  of  operators  on  the  up  and  downgoing  waves  generated  at  the  source. 

Useful  approximations  in  the  operator  development  to  concentrate  attention 
on  the  dominant  modes  of  propagaticn  can  be  made  by  exploiting  analogies 
with  earlier  work  on  stratified  media.  The  particular  implementation  for 
the  operators  can  be  adapted  to  suit  the  particular  nature  of  the  process 
being  described. 

This  operator  representation  has  been  used  to  look  at  the  theoretical 
characteristics  of  a  variety  of  potential  discriminants  between 
earthquakes  and  underground  nuclear  explosions  based  on  the  ratio  of  P  and 
S  wave  amplitudes. 

A  full  description  of  the  operator  approach  and  the  discrimination  results 
are  given  in  the  attached  paper:  Oh  the  nature  of  regional  seismic  phases  ■ 
I  -  phase  r^resen  tat  ions  for  Pn,  Pg,  Sn  and  Lg. 
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Many  existing  techniques  for  describing  the  propagation  of  guided  waves 
in  heterogeneous  nedia  can  be  regarded  as  inplenentations  of  the  operator 
approach.  The  coupled  node  technique  of  Kennett  (1984)  allows  for 
interactions  with  heterogeneity  by  allowing  the  transfer  of  energy  between 
differmt  nodes  with  different  horizontal  wave-nunbers,  thereby 
redistributing  the  energy  between  different  group  velocity  and  so  changing 
the  character  of  the  seisnic  record.  This  approach  is  particularly  suited 
for  exanining  the  effect  of  nodest  heterogeneity  on  regional  S  wavetrains 
which  can  be  described  by  a  United  nunber  of  discrete  nodes. 

In  the  acconpanying  paper;  wave  pn^pagation  in  heterogeneous  media. 
a  sequence  of  n^els  with  varying  levels  of  heterogeneity  have  been  used 
to  detemine  the  nerits  and  linitations  of  the  coupled  n^e  cooputation 
scheae.  The  technique  works  well  with  heterogeneous  nodels  in  which  the 
local  seisnic  velocities  differ  fron  the  stratified  reference  by  up  to 
2  per  cent  without  significant  perturbation  of  the  major  interfaces  (such 
as  the  crust-mantle  boundary).  Localised  changes  in  seismic  velocities 
and  density  can  be  even  larger,  but  it  is  difficult  to  account  for  shifts 
of  more  tt^  2-3  km  in  the  position  of  the  Moho. 

The  nature  of  the  modifications  produced  by  Lg  propagation  through  with 
heterogeneity  can  be  clearly  seen  by  constructing  theoretical  seismograms 
with  a  narrow  frequency  band  and  then  incorporating  the  inter-node 
coupling.  The  shape  of  the  main  amplitude  maximum  is  modified  and  there 
can  be  transfer  of  energy  to  Sn  which  terxis  to  extend  the  duration  of 
significant  signal. 
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On  the  nature  of  regional  seismic  phases  1 
-  phase  representations  for  Pn,  Pg,  Sn,  Lg 

B.L.N.  Kennett 

Research  School  of  Earth  Sciences, 

Australian  National  University,. 

6.P.O.  Bok  4,  Canberra,  A.C.T.  2601 
Australia 

Summary 

I 

An  operator  development  of  the  seismic  wavefield  is  used  to 
generate  descriptions  of  the  propagation  processes  contributing 
to  the  main  regional  seismic  phases  Pn ,  Pg ,  Sn ,  Lg .  These 
operator  forms  are  valid  for  laterally  heterogeneous  crust  and 
mantle  models  and  include  the  major  processes  of  interconversion 
between  wavetypes. 

These  representations  of  the  regional  phases  are  used  to 
examine  the  theoretical  basis  for  discriminants  between 
earthquakes  and  underground  explosions  based  on  the  relative 
amplitudes  of  P  and  S  phases.  The  ratio  of  Sn  to  Pn  amplitude 
looks  promising  as  a  high  frequency  discriminant.  However,  the 
ratio  of  Lg  to  Pn  amplitudes  is  not  as  useful  because  of  the 
complex  nature  of  the  propagation  characteristics  of  Lg. 

Key  Words:  regional  phases,  wave  propagation,  laterally 

heterogeneous  media,  Pn,  Pg ,  Sn,  Lg ,  source  discrimination 


-3- 


Intrcxluction 


The  propagation  of  aeismic  phases  at  regional  distances  has 
been  a  topic  of  continuous  interest  since  the  work  of  Mohorovicic 
(1909)  on  the  Kulpatal  earthguake  in  Croatia,  and  these  phases 
have  been  eKtensively  used  to  determine  the  structure  of  the 
crust  and  uppermost  mantle.  These  structural  studies,  especially 
for  P  waves,  lead  to  the  development  of  many  analysis  tools  such 
as  the  computation  of  theoretical  seismograms  for  realistic 
models  (see  e.g.  Fuchs  *  Muller  1971).  A  major  impetus  to 
understand  the  detailed  characteristics  of  the  regional 
wavefield  has  come  with  efforts  to  monitor  underground  explosions 
of  low  yield  and  discriminate  them  from  earthquakes 
(see  e.g.  Pomeroy  et  al  1962). 

Despite  this  long  term  interest,  there  has  been  considerable 
debate  about  the  nature  of  the  propagation  processes  giving  rise 
to  the  most  prominent  of  the  observed  regional  phases  Pn,  Pg,  Sn, 
Lg.  The  object  of  this  paper  is  to  provide  a  theoretical 
description  of  the  propagation  of  these  regional  phases  which  can 
be  used  as  the  basis  for  discussions  of  the  effects  of  structure 
and  source  type.  The  development  is  based  on  the  operator  methods 
of  Kennett  (1984),  to  get  a  representation  of  the  full  wavefield 
in  crust  and  mantle  models  with  lateral  variation  in  seismic 
properties.  The  operator  results  are  then  broken  down  into 
contributions  associated  with  particular  wavetypes  using 
approximations  introduced  by  Kennett  (1986a).  The  net  result  is  a 
specific  description  of  the  propagation  processes  contributing  to 
Pg,  Lg  and  Pn ,  Sn  with  inclusion  of  the  most  important 
interconversions  between  wavetypes,  for  a  laterally  heterogeneous 
struc  ture . 

These  representations  of  the  regional  phases  are  tnen  used  to 
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examine  the  behaviour  of  some  proposed  discriminants  between 
earthquakes  and  underground  explosions  based  on  the  relative 
amplitudes  of  P  and  S  phases.  The  preferred  phase  for  use  for  P 
waves  is  t'le  Pn  phase  because  it  represents  a  distinct  arrival 
at  the  beginning  of  the  record.  The  size  of  the  later  Pg  phase  is 
too  dependent  on  the  nature  of  the  propagation  path  to  be  general 
useful.  The  ratio  of  the  amplitudes  of  the  Pn  and  Sn  phases 

looks  promising  as  a  high  frequency  discriminant  for  events 

/ 

which  are  large  enough  fon  Sn  to  be  observed 

clearly.  However,  the  ratio  of  the  amplitudes  of  Lg  and  Pn 

is  not  as  useful  because  of  the  complex  nature  of  the 

of  the  propagation  characteristics  of  Lg  in  heterogeneous  media. 

Computations  based  on  an' implementation  of  the  operator 
representation  for  laterally  varying  madia  via  wavenumber  mixing 
in  the  transform  domain  will  be  presented  in  a  companion  paper. 

2.  General  description  of  regional  phases 

Me  will  here  show  how  the  operator  development  of  the  seismic 
wavefield  introduced  by  Kennett  (1984,  1986a)  can  be  adapted  to 
the  description  of  regional  seismic  phases. 

2.1  Reflection  end  trenmmimeiiin  operators 

Consider  a  heterogeneous  region  bounded  by  smooth  surfaces  A 
(  z  ■  f.  fx  J,  X  »  )  and  C  (  z  »  f-  fx,)  )  which  do  not 

deviate  far  from  horizontal  planes.  In  oriirr  to  isolate  the 
region  from  its  surroundings,  we  consider  it  to  be  bordered  on 
each  side  by  uniform  media  a,  c  (see  fig  1).  Then,  it  there  is 
a  downgoing  wavefield  D*  incident  from  the  uniform  medium  a 
upon  the  region  AC,  the  reflected  field  back  into  a  is  written. 
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a>  *  > 


(2.1) 


Mhcr*  R 


i*  «  rv-flvction  op«r«tor  for  downgoing  waves;  and 


the  transmitted  field  in  c  is  represented  as 


U  (M  )  m  iO*  >  . 


(2.2) 


4e 


in  terms  of  a  transmission  operator  T  _  R  _  ,  T 

^  J) 

include  the  effect  of  any  multiple  reflections  within  the  region 
AC.  For  an  incident  upgoing  wave  from  the  uniform  medium  c  we 
may  similarly  introduce  the  corresponding  operators  ^  v* 

for  upgoing  waves. 

These  reflection  and  transmission  operators  are  not  tied  to  any 
particular  computational  implementation,  and  include  the  position 
dependence  of  the  reflected  and  transmitted  wavefields.  In  a 
stratified  medium  the  action  of  R^  ,  can  be  evaluated 

via  matrix  multiplications  in  the  frequency-wavenumber  domain  and 
subsequent  inversion  of  the  transforms.  More  complex  reflection 
and  transmission  processes  can  be  represented  by  including  the 
action  of  further  operators.  The  operator  algebra  is 
noncommutative  and  operators  always  act  to  their  right. 

Now  divide  the  region  AC  by  the  introduction  of  an  intermediate 
surface  B  )  >  into  two  subregions  AB  and  BC.  The 

composite  reflection  and  transmission  operators  R  ,  T'*' 
can  be  represented  in  terms  of  those  for  the  subregions  as 


1^*  ♦  t"**  r  I  -  r'*"  R*‘  J-’  T**  . 

*  "a  v»  » 


T*'  r  1  -  r'’*  R*'  r’  t** 

D  u  a  a 


(2.3) 


where  I  is  the  identity  operator.  We  recall  that  operators  act 
to  their  right  and  so  the  physical  significance  of  each 
expression  can  be  obtained  by  reading  from  right  to  left. 
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The  operator  inverse  t  I  ~  l”'  represents  ell 

ir ' erections  between  the  regions  AB  end  BC  end  hence  ell 

AS  $C  • 

multiple  reflections  within  AC.  We  will  refer  to  [  I  -  R^  3* 

es  the  reverberetion  operetor  for  the  region  AC. 

The  operetor  elgebre  hes  the  seme  structural  form  es  the 
reflection  metrix  methods  foe  stratified  media  (see  e.g.  Kennett 
1963).  This  correspondence  arises  because  the  reflection  matrices 
in  the  transform  domain  ere  just  one  implementation  of  the 
operator  forms.  As  a  result,  we  ere  able  to  exploit  e  wide  range 
of  results  end  approximetions  from  the  stratified  case  with  a 
reinterpretation  of  matrices  as  operators. 

We  have  so  far  considered  the  reflection  end  transmission 
operators  as  acting  on  the.  entire  incident  wavefield.  But  within 
the  uniform  media  we  have  supposed  to  border  the  region  of 
interest,  we  can  split  up  the  wave  field  into  three  independent 
wave  types.  In  isotropic  media  these  will  be  P-waves,  SV-waves 
end  SH-waves,  but  in  anisotropic  media  the  decomposition  by 
wevetype  may  not  lead  to  readily  identifiable  physical  character. 
We  will  however  label  the  three  wave  types  by  P,  S  and  H,  as  in 
Kennett  (I9a6a)  and  allow  for  the  interconversions  between 
wavetypes  arising  from  the  nature  of  the  medium.  We  can  partition 
the  operators  by  input  and  output  wavetypes  so  that  we  can  write 


e.g. 


R 


AC 

s 


(2.4) 


where  we  have  adapted  the  convention  that  R represents 
scattering  from  wavetype  2  into  wavetype  1  (Kennett  1983).  With 
this  convention  the  partitioned  forms  of  composite  operators  can 
be  constructed  by  equivalent  rules  to  matrix  multiplication,  for 
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CKample  the  PP  element  ofR  R  is  given  by 


(r""  R^)  -  r 

'  w  »  life 


dc  /4a  8c 

^  R  R 
ofe  see  ves  *»>• 


R  R 

Ufa  »He 


(2.5) 


Even  for  en  isotropic  medium,  the  action  of  heterogeneity  will 
often  be  to  produce  out  of  plane  scattering  and  so  couple  the 
P-SV  and  SH  fields  together. 


2,2  Opmrmtor  dmmcriptian  of  thm  seismic  leavefieid 

In  the  regional  phase  situation  we  consider  a  source  lying  in 
the  crustal  2one  at  a  level  S  Iz  *  )  with  the  crust/mantle 

boundary  defined  by  the  surface  C  ( z .  ^  fig  2. 

If  the  source  was  placed  in  a  uniform  medium,  with  the  local 
properties  at  the  source,  it  would  radiate  upward  and  downward 
contributions  U*  ,  D*  .  By  splitting  the  medium  at  the 
source  level,  we  can  work  in  terms  of  reflection  and  transmission 
operators  above  and  below  the  source  level  S  and  introduce  the 
source  via  its  upward  iU*  )  and  downward  ID*  )  radiation  terms. 

The  resulting  representation  of  the  seismic  wavefield  takes  the 
form 

U  ;  -  fi  T^*  C  1  -  «**■  R^*  J'*  (O*  *  R**"  P*  ),  (2.6) 

where  SL  indicates  the  entire  region  below  the  source  level  and  fS 

indicates  the  zone  above  the  source  including  the  free  surface. 

represents  the  transmission  of  an  upgoing  wave  from  the 

source  level  up  to  the  free  surface  and  H  is  the  surface 

amplification  factor  arising  from  the  interference  of  up  and 

downgoing  wavefields  at  the  free  surface.  The  inverse  operator 
SL  fs  -1 

Cl-  R^  R^  ]  includes  all  multiple  interactions  between 
the  regions  above  and  below  the  source.  Equation  (2.6)  for  the 
seismic  wavefield  can  be  rewritten  in  a  way  that  emphasises  the 
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surface  reflected  phases 


u  (X  ) 

/V  ^  ^  ^ 


w 

l  uU 

fs 

*  M 

T 

u 

C  I  - 

(2.7) 


In  this  representation  for  the  wavefield  we  have  the  reverberation 
operator  for  the  entire  structure  combined  with  reflection  from 
beneath  the  source  lev^l  .  Since  we  have  the  operator  identity 


C  I 


a*-  ft  .1  at 

R..  1  R. 


^Si.  fS  CL  .4 


we  can  view  the  reverberation  sequence  as  occurring  on  either  the 

.St 


source  or  the  receiver  sides  of  the  main  reflection 


Me  can  now  begin  to  separate  phases  propagating  in  the 


crust  from  those  with  a  mantle  component  by  splitting  the 


reflection  operator  R  just  below  the  crust/mantle  boundary  C 


using  (2.3) 
S'- 


R**  *  T**  [  I 

i  u 


_ct  _ae  _ci.  _*e 

R  R  1  R 


(2.8) 


This  separation  includes  reflection  of  downward  travelling  waves 

.ac 


at  the  crust-mantle  interface  in  the  crustal  operator  .  The 

dominant  multiples  in  regional  propagation  arise  from  the  free 
surface,  but  there  can  be  a  contribution  from  a  complex 
crust-mantle  transition.  We  will  therefore  rewrite  (2.8)  in  a 
form  which  emphasises  the  separation  into  crust  and  mantle 
operators : 


^  m.  SC 

T  R  T  ^ 
u  ^ 


(2.9) 


Here  R^  includes  all  mantle  reflection  effects,  both  direct 


reflection  (R^  )  and  any  modulation  from  the  effect  of  the 


CL  sc  ~ 

structure  of  the  crust/mantle  transition  through  t^“RjjRy  3 


We  would  like  to  make  use  of  the  expansion  (2.9)  to  separate 


i 
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multiple  reflections  in  the  crustal  zone  from  the  overall 

*le  make  use  of 

the  operator  identity 


SL  fs  , 

reverberation  operator  t  I  ~  ^  ^  J 


C  I  -  A  -  B  ] 


CI-A3*^  ♦CI-A]'"bCI-A-B]‘ 


which  leads  to  an  espansion  In  multiple  powers  of  the  operator 
combination  B  [  1  —  A  3*^  .  With  the  identification  of  A  with 
R*  R^*  and  the  mantle  contribution  T*  R.  R^*  as  B,  we 

can  develop  an  eKpansion  df  the  form 


SL  fi  -1 

C  1  -  R,  «  „  ] 


c « -  «r  3  * 


♦  I  1  -  R^  R^*  3-'  t“  r'^ 


T*R^*t  I  “  R^  A^*  3  * 

(2.10) 


and  we  can  recognise  t  1  -  R^  R^*  3'*  *»  the  reverberation 
operator  for  the  crust. 

When  we  make  use  of  the  crust~mantle  decomposition  in  (2.9)  and 
(2.10)  in  the  expansion  for  the  seismic  wave  field  (2.7)  we  can 
set  up  a  hierachy  of  propagation  terms  emphasising  the  different 
classes  of  wave  interaction  with  the  structure  of  significance  for 
regional  phases 


a,  ^ 


w  T^*  iu*  1 


W  T^*  Cl-  R**  R^*  3'*R*CB*  ♦  R^  1/ *  3 


(2.11a) 


(2.11b) 


..  S  T^*  r  .  T*  T*  b’ 

♦  W  ^  ~  'v  "j 


.to*  *  R^*  o*  ) 

•  V  • 


t  I  -  R^  r"  3-‘ 


(2.11c) 


The  contribution  (2.11a)  represents  the  direct  waves  transmitted 
upward  from  the  source  level.  The  part  (2.11b)  represents  the 
crustal ly  trapped  waves  such  as  Pg.  Lg  involving  multiple 
interactions  with  the  free  surface  and  the  crust/mantle  interface. 
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The  contribution  (2.11c)  represents  dominantly  upper  mantle 
propagation  such  as  Pn,Sn  with  the  possibility  of  neai — source  or 
near-receiver  crustal  reverberations.  The  higher  order  terms 
involve  multiple  mantle  interaction. 

Equation  (2.11)  thus  provides  an  operational  description  of 
the  main  regional  seismic  phases  in  terms  of  the  propagation 
characteristics  of  the  regions  above  and  below  the  source.  For 
shallow  sources  the  natural  divisions  of  the  heterogeneity 
structure  may  well  not  coincide  with  the  placing  of  the  source. 
However,  with  a  little  extra  effort  we  can  generate  an  operator 
description  for  fixed  upper  and  lower  crustal  zonesi  the  details 
are  presented  in  the  Appendix.  For  heterogeneity  superimposed  on 
stratification  the  operators  can  be  simulated  by  wavenumber 
coupling  techniques  in  the  frequency-wavenumber  domain  (Kennett 
1986).  Equation  (2.11)  can  then  be  used  to  synthesise  the 
response  for  the  various  phases  including,  for  example,  different 
scales  of  heterogeneity  in  the  mantle  and  the  parts  of  the  crust 
lying  above  and  below  the  source. 

However,  (2.11)  provides  a  composite  description  of  the  entire 
seismic  wavefield  and  provides  no  immediate  separation  into  the  P 
and  S  regional  phases.  In  order  to  make  a  specific  description  of 
the  individual  phases  Pg,  Pn,  Lg,  Sn  we  need  to  look  in  detail  at 
near  source  and  near  receiver  processes  and  the  character  of  the 
crustal  reverberations. 

2.1  Meat — source  effects 

The  action  of  the  source  enters  the  operator  description  of  the 

<  /s  , 

wavefield  through  iD  +  U  >  which  represents  the  combination 

of  the  direct  downward  radiation  from  the  source  and  reflections 
from  above  the  source,  which  will  be  dominated  by  the  effect  of 
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th«  frm»  «urf«cc. 

In  t»rms  of  the  thro*  b«ftic  M«v»typ*«  (P,S|H)  we  have  the 
eKplicit  e»pre««lon« 


4- 
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4- 
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U*  - 

Ot 

4* 

< 

4> 

a<« 

", 

4- 

"x  X 
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MS 

, 

Mil 

••  ' 

(2.12) 

allowing  for  full  interconversion  above  the  source. 

/ 

For  an  OKplosive  souhcei  the  doeinant  radiation  will  be  as  P 
in  ,  Up  and  S  wave  contributions  will  be  generated  from  free 
surface  reflections!  these  would  normally  be  mainly  6V,  but  some 
SH  waves  could  occur  due  to  heterogeneity  above  the  source.  For  a 
large  explosion,  there  is  also  the  possibility  of  tectonic 
release  close  to  the  source  inducing  direct  S  radiation.  Similar 
effects  will  arise  from  extended  or  multiple  sources.  Thus  if  we 
consider  the  SV  wave  radiation  from  an  explosion  we  can  group 
the  source  contributions  in  terms  of  their  likely  significance! 


MX 


U„ 


(2.13) 


the  first  term  would  be  generated  with  a  stratified  medium  code 
and  a  simple  source  model,  the  second  term  is  associated  with 
tectonic  release  and  the  third  would  arise  almost  entirely  from 
heterogeneity  effects. 

For  an  earthquake  source,  we  have  to  take  account  of  P,SV  and 
SH  radiation  at  the  source  itself  and  heterogeneity  will  have  the 
effect  of  linking  the  SV  and  SH  wavefields.  In  this  context  we 
should  note  that  such  coupling  would  be  predicted  directly  for 
anisotropic  models.  However  if  we  adopt  an  isotropic  reference 
model,  the  action  of  anisotropy  and  heterogeneity  are  formally 
equivalent  and  combined  in  the  same  operator  description  (although 
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the  nature  of  the  coupling  terms  will  be  different). 

As  shown  in  the  Appendix,  the  effect  of  heterogeneity  in  the 
source  zone  will  be  to  modify  the  effective  source  to 
X*  U*  },  where  X*  represents  multiple  scattering 

effects  in  the  source  region.  The  action  of  X*  will  be  to 
introduce  further  coupling  between  wavetypes.  This  will  have  the 
effect,  for  example  of  generating  apparent  downward  radiation  of  S 
waves  from  even  a  simple  explosion  due  to  heterogeneity  induced 
contributions  of  the  type  'Xf  D  ,  xf  Z>  . 


2.2  NsMt — recaiv^r  effects 

The  receivers  are  normally  situated  at  or  near  the  free  surface 

and  so  lie  in  a  region  of  concentrated  small-scale  heterogeneity 

even  in  regions  with  apparently  homogeneous  geology.  Signal 

generated  noise  from  the  vicinity  of  the  receivers  is  frequently 

one  of  the  major  contaminants  of  high  frequency  seismic  data. 

In  the  operational  description  of  the  seismic  phases  (2.11) 

"  fs 

such  effects  are  included  in  the  term  W  T^  .  For  low 
frequency  waves  this  operator  will  be  close  to  that  for  a 
stratified  medium.  But,  at  high  frequencies  we  will  have  to  allow 
for  significant  deviations  from  the  stratified  results  induced  by 
three-dimensional  velocity  structure  and  local  topography  near  the 
receivers.  We  will  therefore  need  to  assume  that  any  incident 
wavetype  at  the  level  S  can  produce  all  three  components  of 
motion  at  the  surface,  and  that  wavetype  coupling  is  likely  to  be 
significant  (both  in-plane  e.g.  P-SV  and  out-of-plane  e.g.  P-SH). 
Thus  we  will  write 


i 

Wzs 

Wzx  \ 

= 

i  W 

M 

f 

w 

V  ^ 

r5 

TM  / 
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in  t»rm*  of  displacement  components  on  the  vertical  (Z)>  radial  (R) 
and  tangential  (T)  components  at  a  receiver,  for  a  particular 
wavetype  incident  at  the  level  S.  The  character  of  the  surface 
motion  may  eell  include  other  wavetypes  than  the  incident.  For  a 
horizontally  stratified  medium  are  aluays  zero  and 

would  only  be  significant  in  the  near-field  of  a  source. 


2.3  Rmvmrbmr^tian  mKp^mx^n 

As  we  have  seen  reverberation  operators  of  the  type 
%€  4S  mi 

C  I  ~  R^  ]  play  an  important  role  in  the  description  of 
the  seismic  wavefield  at  regional  ranges.  The  compact 
representation  of  the  operator  inverse  does  not  allow  the 
identification  of  reverberation  sequences  associated  with 
individual  wavetypes. 

The  crustal  reverberation  operator  is  the  inverse  of 
Cl-R^Ry  ]  which  can  be  written  in  partitioned  form  as 


t  I  _  J 


Xn.  -  V  - 


-  Y. 


-  Y. 


\  -  Vw 

where  the  diagonal  terms  are  e.g. 

fs 


(2.15) 


Y  X  -  Y  / 

H*  »«  HM  / 


R**  r" 

art  Mr 


n  •*<=  r, 

R  R 

art  UMP 


(2.16) 


and  the  terms  involving  conversions  take  the  form 


Y 


rt 


♦ 


R 


fs 


Using  (2.15),  it  is  possible  to  construct  eicact  expressions 
the  partitions  of  the  reverberation  operator  in  terms  of  the 
individual  wavetypes  (Kennett  1786,  section  5.2).  However, 


(2.17) 

for 
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the  resulting  expressions  are  rather  complex,  especially  Mhen 
interaction  between  all  three  wavetypes  are  included. 
Fortunately,  multiple  conversions  between  wavetypes  are  unlikely 
to  be  of  very  great  signif icance,  and  so  we  may  employ  a 
convenient  approximation  to  the  reverberation  operator  in  which 
no  more  than  two  conversions , are  considered! 


C  I 
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(2.18) 


In  this  approximation  we  have  been  able  to  extract  the 

-4  ^.4  -4 

reverberation  terms  for  the  individual  wavetypes  ,  X^  ,  X^,^  . 
Conversions  between  wavetypes  are  accompanied  by  the  reverberation 
sequences  for  both  primary  and  secondary  wavetypes. 


3.  Representation  of  specific  seiseic  phases 

In  the  previous  section  we  have  seen  how  we  may  describe 
the  seismic  wavefield  at  regional  ranges  with  an  operational 
development  and  have  also  considered  near-source,  near-receiver 
and  crustal  reverberation  effects.  With  these  preliminaries,  we 
are  in  a  position  to  produce  representations  of  specific  seismic 
phases.  We  will  start  with  crustal ly  guided  waves  and  then  turn 
attention  to  those  phases  whose  paths  lie  mainly  in  the  mantle. 

J.i  Crustal ly  guided  waves 

The  first  group  of  waves  whose  paths  lie  in  the  crust  is  that 
associated  with  transmission  upward  from  the  source  (2.11a) 

u'*=WT^‘{t/’)  (3.1) 

—  •  u  .. 

These  direct  P  or  S  waves  are  important  close  to  the  source,  but 
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at  greater  distances  are  reduced  in  amplitude  by  geometrical 
spreading.  Since  the  direct  waves  travel  through  the  lower 
velocity  material  near  the  surface,  they  are  overtaken  at  larger 
distances  by  phases  travelling  mostly  through  the  higher  velocity 
material  at  depth,  and  then  play  a  fairly  minor  role  in  the  coda. 
Of  much  greater  significance  are  those  waves  which  have  been 
reflected  back  from  beneath  the  source  level  and  than  have  been 
involved  with  reverberation  through  the  whole  crustal  zone  (2.11b) 


C  I  -  R*  ♦  R^*  W*  > 


We  recall  that  the  level  C  lies  just  below  the  crust-mantle 

sc 

transition  so  that  R^  includes  reflection  from  this  interface 
(e.g.  PmP,  SmS  arrivals).  A  schematic  representation  of  the 
various  operator  interactions  in  (3.2)  is  shown  in  figure  3. 

At  low  frequencies,  the  structure  beneath  the  interface  will  be 
involved  in  the  reflection  process  and  so  influence  the  crustal 
reverberations  (Burdick  &  Helmberger  1988).  However,  for 
frequencies  higher  than  about  0.5  Hz,  the  representation  (3.2) 
will  include  the  dominant  contribution  to  the  crustal ly  guided 
phases  Pg,(.g.  The  effect  of  source  depth  will  enter  into  the 
CMpression  (3.2)  through  the  compound  source  term  iD*  *  R^'  U*  }• 
The  balance  between  downward  and  upward  waves  will  be  of 
particular  relevance  for  directional  sources  such  as  earthquakes. 

With  the  aid  of  the  decomposition  of  the  crustal  reverberation 
operator  by  wavetypes  discussed  in  section  2.3,  we  can  now  look 
at  the  individual  seismic  phases. 


3.J.i  Pg 

The  Pg  phase  is  composed  primarily  of  P  waves  guided  in 
the  crust  although  the  coda  is  likely  to  have  an  increasing 
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proportion  of  converted  wave&.  Conversion  near  the  source  will 
normally  be  significant  in  transferring  S  energy  frc.  the  source 
radiation  to  the  P  wavefield.  The  base  of  the  weathered  zone  near 
the  surface  will  be  one  of  the  most  likely  places  for  conversion 
from  P  to  other  wavetypes  and  so  transfer  energy  to  the  tangential 
component.  These  effects  can  be  included  by  making  a  full 
representation  of  the  transmission  through  the  surface  zone  (2.14). 
The  most  significant  contribution  to  the  Pg  phase  will  be  from 
those  wave  propagation  prgcesses  with  just  P  reverberations  in  the 
crust  With  the  form 


C  R 
+R 
+R 


♦ 

i/j 

(o. 

+ 

c: 

MC 

♦ 

)1 

*  > 
) 


(3.3) 


where  is  the  crustal  reverberation  operator  for  P  waves 


r  1  -  o*  o" 


(3.4) 


is  the  transmission  operator  for  the  surface  motion  generated 
by  an  incident  upward  P  wave  at  the  source  level 


M  ■  -  C  ]  , 


(3.5) 


this  allows  for  scattering  of  energy  out-of-plane  and  possible 
conversions  near  the  surface. 

The  contributions  to  the  coda  of  Pg  involve  multiple 
reverberation  trains  such  as 

x'’y  x"’  =  Cl-  R**  R^‘  ]*’  (R*'  R^*  +R'*'  )tl-  Rf_R^t»^’*. 

fr  rr  FT  imr  »w.  v»r  »*«  vhp  • 

(3.6) 

and  similar  terms  with  conversion  between  reverberations  of  the 
type 


Cl-  R**"  r”  }■’ 

DTP  ppp 


(3.7) 
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with  thorn*  proc**s  involving  SH  w«v*s  lik*ly  to  be  of  least 
•ignif icanc*<  Sine*  all  th*s*  contributions  require  conversions 
between  wavetypest  they  are  likely  to  be  individually  small  but 
the  cumulative  effect  of  a  nvtmber  of  such  terms  can  still  be 
noticeable. 

The  character  of  the  Pg  wavetrain  will  be  strongly  dependent  on 
the  nature  of  the  crustal , reverberation  operator  Multiple 

reflections  in  the  crust  for  P  provide  only  partial  trapping  of 
energy,  because  of  conversion  from  P  to  6  at  the  surface  and  the 
possibility  of  transmission  through  the  crust-mantle  transition 
into  S  waves  in  the  mantle. 

For  laterally  hoaiogeneous  crustal  models,  the  detailed  studies 
of  Olsen  et  al  (1983),  Campillo  *t  al  (1984)  have  shown  how  the 
Pg  phase  is  built  up  from  multiple  PmP  reflections.  However,  it 
is  difficult  to  sustain  the  Pg  phase  to  even  500  km  range  if 
there  is  not  low  velocity  material  near  the  surface.  A  crustal 
low  velocity  zone  may  help  to  maintain  Pg  amplitudes  by 
interference  of  multiples  separated  by  only  a  short  time  delay. 

In  a  shield  type  model,  the  Pg  energy  leaks  rapidly  into  5 
with  each  surface  reflection  and  the  Pg  phase  tends  not  to  be 
observed  at  large  ranges.  The  presence  of  lateral  heterogeneity  will 
tend  to  enhance  the  decay  of  Pg  with  range,  since  scattering  will 
generally  remove  energy  from  the  waveguide.  These  results  help  to 
explain  why  Pg  is  more  readily  observed  in  the  western  United 
States,  where  surface  P  wave  velocities  tend  to  be  low,  than  in 
the  eastern  United  States  where  older  geological  structures  bring 
higher  vel  ocity  material  close  to  the  surface. 


-18- 


17 


J.1.2  Sg,  Lg 

At  close  ranges,  the  crustal ly  guided  S  waves,  analogous  to  Pg , 
are  commonly  referred  to  as  Sg.  But  at  greater  distances  as  the 
wavetrains  become  more  complex,  the  notation  Lg  is  frequently 
used.  Originally  the  term  Lg  was  applied  to  intermediate 
frequency  arrivals  with  a  group  velocity  around  3.5  km/s  (Press 
&  Ewing  1952),  but  more  recently  the  usage  has  been  transferred  to 
,  higher  frequency  waves  (  >  1  Hz). 

When  considering  these  ,S  wavetrains  we  have  to  take  account  of 
the  different  classes  of  wave  propagation  processes  associated 
with  the  two  S  wavetypes.  In  a  horizontally  stratified, 
isotropic,  medium  the  SV  waves  would  propagate  in  a  vertical  plane 
and  SH  in  a  horizontal  plane.  The  SV  waves  part  of  Lg  can  then  be 
described  by  a  superposition  of  many  higher  modes  of  Rayleigh 
waves  and  the  SH  wave  part  by  a  similar  sum  of  Love  wave  terms. 
However,  in  a  real  medium  there  will  be  cross-coupling  between  the 
wavetypes  due  to  out-of-plane  scattering  arising  from 
heterogeneity  in  the  crustal  waveguide. 

The  phase  velocities  of  importance  for  guided  S  waves  are  such 
that  P  waves  will  only  propagate  in  the  near  surface  region  (if  at 
all).  As  a  result,  we  only  need  to  retain  P  waves  in  the 
description  of  near — source  and  near — receiver  processes.  Within  the 
full  crustal  waveguide  the  behaviour  will  be  dominated  by  S  wave 
reverberations,  but  for  each  wavetype  we  must  make  a  full 
allowance  for  near-receiver  effects  including  coupling  to  P  waves 
and  Rg  waves  (fundamental  mode  Rayleigh  waves)  largely  confined  to 
the  near  surface. 

In  order  to  simplify  the  notation  we  introduce 
representing  the  net  downward  radiation  for  SV ,  SH  waves 
including  the  effect  of  reflection  processes  above  the  source 
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Th«  demin«nt  rmtlmctLon  proc««Ms  «bov*  the  source  are  likely  to 
be  those  without  change  of  wavetype.  The  main  effects  of  the 
guiding  process  will  be  described  by  the  crustal  reverberation 
operators  for  the  two  S  wave  types  representing  the  effects  of 
multiple  internal  reflections  within  the  crust.  Using  the 
notation  of  section  2.3  we  will  write 
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Within  the  crustal  waveguide  we  can  expect  cross-coupling  between 
wavetypes  through  the  terms  Y^j,  ,  Y^^  ,  Y^  with  principal 
contributions 
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and  terms  of  lower  importance  involving  P  wave  conversion.  We 
must  allow  for  full  coupling  in  propagation  through  the 
near-surface  xone.  Tor  incident  S  waves  the  two  transmission 
terms  can  be  represented  as 
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on  the  nature  of  the  main  sequence  of  multiple  5  wave  reflections. 
Within  the  waveguide,  cross-coupling  between  wavetypes  will 
increase  the  complexity  of  the  wavepropaga tion  process.  With 
dominantly  SV  reverberat ion ,  the  major  contributions  to  Lq  can  be 
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with  a  similar  form  for  the  dominantly  SH  wave  propagation 
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Because  the  wavespeeds  of  the  different  wavetypes  are  so  similar, 
we  need  to  retain  a  more  complete  representation  than  for  Pg. 
However,  the  main  effects  will  arise  from  terms  with  a  single 
reverberation  sequence. 

For  the  Lg  phases,  the  trapping  of  S  energy  is  quite  efficient. 
In  the  laterally  homogeneous  case,  there  will  be  total  reflection 
of  S  waves  at  the  surface  and  at  a  crust-mantle  interface;  some 
slight  energy  loss  may  arise  at  the  base  of  the  crust  for  a 
transitional  structure  into  the  mantle.  Even  in  the  presence  of 
lateral  heterogeneity  most  of  the  energy  will  be  trapped.  There 
will  be  some  energy  loss  by  anclastic  attenuation  and  scattering 
and  so  the  reverberation  sequences  ,  X’*^  will  decay  in 
time.  This  leads  to  a  general  reduction  of  the  coda  amplitude 
along  the  wavetrain  as  time  progresses,  whose  envelope  can  usually 
be  appro¥imated  by  an  exponential  to  define  a  coda  Q  (Herrmann 
1980) . 

There  are  a  number  of  regions  where  the  propagation  of  the  Lg 
wave  IS  interrupted  by  some  tectonic  feature  e.g.  the  graben 
structure  in  the  North  Sea  (Gregersen  1904)  and  the  mountain  belts 
of  central  Asia  (Ruzaikin  et  al,  1977).  This  disruption  of  the 
wavetrain  is  generally  associated  with  upsetting  the  constructive 
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intarf*r«nc*  of  the  multiple  crustal  S  reflections  comprising 
X jj  ,  xj^  .  As  illustrated  by  Kennett  (19B6b)i  modifications 
of  the  shape  of  the  crustal  waveguide  tend  to  break  up  the 
constructive  interference  and  these  trends  will  be  enhanced 
by  significant  horizontal  velocity  gradients-  In  addition 
anoealously  low  Q  may  be  needed  to  finally  extinguish  the  guided 
energy  (Maupin  i989). 


3.2  Hantim 

Once  we  get  beyond  the' critical  distance  for  reflection  from 
the  crust-mantle  transition,  we  have  the  possibility  of  phases 
with  a  significant  path  in  the  mantle  associated  with  operator 
(2.9)  ,  which  includes  all  reflection  from  beneath  our 
separation  level  C  (just  below  the  crust-mantle  inter  ce). 
The  general  description  of  such  mantle  phases  is  provided  by 
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which  includes  the  possibility  of  reverberations  in  the  crust 
near  the  source  or  near  the  receiver.  We  recall  that  R.  is  a 
contracted  notation 


designed  to  include  all  mantle  reflections  including  those 
arising  from  the  structure  of  the  crust-mantle  transition  zone. 

As  in  (3.2)  the  effect  of  source  depth  enters  into  the  expression 
for  the  wavefield  through  the  compound  source  term  (0‘  +  R^' ^  5* 
We  can  extract  the  main  term  corresponding  to  Pn ,  Sn  from 
(3.13)  as 
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which  includes  the  possibility  of  reflections  from  the  surface  at 
the  source  e.g.  pPn ,  sSn  (see  fig  4).  The  coda  of  these  two 
phases  will  include  crustal  multiples  such  as  PPn,  SSn. 


J.  2. 1  Pn 

The  beginning  of  the  Pn  wave.train  will  be  dominated  by  P  waves 
radiated  downward  from  the  source,  with  some  reinforcement 
slightly  later  from  surface  reflected  phases  (including  S  to  P 
conversion)  for  shallow  sources.  The  mantle  leg  will  be  as  a  P 
wave  so  that  we  may  represent  the  wavetrain  as 
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(3.15) 


where  cross-coupling  between  wavetypes  in  the  near — surface  2one  is 
included  via  the  transmission  term  (3.5). 

Rt  ranges  beyond  300-400  km,  crustal  multiples  of  Pn  begin  to 
contribute  to  the  coda  of  Pn  as  indicated  in  the  representation 
(3.13).  For  the  range  of  phase  velocites  of  importance  for  this 
phase  8. 3-7. 6  km/s,  conversion  of  P  to  5  waves  at  the  surface  is 
more  efficient  than  reflection  of  P  to  itself  so  that  some  crustal 
legs  may  be  as  S  waves.  As  a  result  the  Pn  coda  will  begin  to 
acquire  a  partial  S  wave  character. 


3.2.2  Sn 

The  principal  phase  velocity  interval  for  Sn  (4. 8-4. 3  km/s) 
will  involve  only  limited  interaction  with  P  waves  in  the  near 
surface  zone  and  so  we  only  need  to  consider  S  wave  propagation 
over  the  bulk  of  the  path.  The  most  significant  other  class  of 
wavetype  conversion  is  likely  to  be  between  the  two  S  waves, 
induced  by  the  presence  of  heterogeneity.  We  will  therefore 
compound  the  transmission  up  and  down  below  the  source  level  with 
the  mantle  reflection  processes  into  a  composite  operator 
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and  than  partition  tha  whola  procas*  by  incidant  and  amargant 
Mavatypa.  Than,  with  tha  aid  of  tha  compact  effactiva  sourca 
notation  introducad  in  saction  3.1.2  we  can  raprasant  the 
dominantly  SV  wava  contribution  to  8n  as 
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with  an  aquivalant  form  for  tha  SH  wave  contribution 
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As  for  Pn,  tha  coda  of  Sn  will  includa  crustal  multiples.  Since 
reflection  at  tha  surface  involves  only  limited  coupling  to  P, 
this  will  be  a  relatively  efficient  process  and  as  such  a  major 
contributor  to  the  character  of  Sn  at  larger  ranges  (Kennett 
1985). 

Although  the  various  operator  descriptions  take  similar  forms 
for  tha  crustal ly  guided  and  mantle  waves,  we  must  recall  that  the 
ranges  of  phase  velocity  for  which  the  representations  are  useful 
are  rather  different.  We  have  been  able  to  partition  the  seismic 
wavefield  via  the  dominant  modes  of  propagation,  but  we  must 
recall  that  the  phases  do  not  exist  in  isolation.  Thus  the  Sn  and 
Lg  arrivals  are  superimposed  on  the  Pg  coda  and  the  spectral 
character  of  the  actual  wavefield  can  be  dominated  by  coda 
contamination  (Blandford  1980).  Also,  we  have  included  all  mantle 
propagation  effects  within  the  Pn,Sn  representations,  energy 
return  from  structure  well  below  the  crust/mantle  interface  can  be 
of  considerable  importance  at  larger  ranges  and  may  well  occur  in 
a  time  window  that,  for  an  isolated  record,  could  be  confused  with 
the  crustal  multiples  of  shallower  propagating  energy. 
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4.  P/S  amplitude  ratios  for  regional  phases 

With  these  relatively  compact  descriptions  of  the  regional 
seismic  phases,  we  are  now  in  a  position  to  look  at  the  behaviour 
of  possible  discriminants  between  different  types  of  seismic 
sources.  For  discrimination  between  earthquakes  and  enplosions, 
such  criteria  are  based  on  the  observation  that  earthquakes  should 
be  much  more  efficient  generators  of  S  waves  than  explosions,  as 
demonstrated  very  elegantly  by  Gilbert  (1973).  For  regional 
phases  this  result  implies  that  we  would,  in  general,  expect 
larger  amplitude  Sn  and  Lg  phases  relative  to  the  size  of  the  P 
waves  for  earthquakes  than  explosions.  Some  earthquake 
observations  near  nodes  in  the  source  radiation  patterns  could  be 
deceptive,  but  with  a  reasonable  coverage  of  azimuths  such 
problems  should  be  avoided. 

Such  discrimination  criteria  based  on  the  relative  amplitudes 
of  P  and  S  regional  phases  have  been  investigated  for  some  time 
but  their  performance  has  been  somewhat  mixed.  Blandford  (1980) 
argues  strongly  for  the  the  use  of  the  ratio  of  the  maximum 
amplitude  before  Sn  (Pmax)  to  the  maximum  after  Sn  (Lg), 
particularly  at  high  frequencies.  However,  Pomeroy  et  al  (1982) 
report  on  a  range  of  other  observations  which  are  not  as 
encouraging.  Clearly,  it  is  desirable  to  look  at  the  theoretical 
basis  for  the  discriminant  to  see  if  we  can  understand  where  the 
problems  arise. 

4.1  Sn/Pn  amplitude  ratio 

This  amplitude  ratio  is  designed  to  compare  the  Pn  and  Sn 
phases  which  both  have  a  significant  mantle  component  to  the 
propagation  path.  Pn  has  the  advantage  that  it  appears  as  a  first 
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Arrival  «nd  only  contonds  with  th*  aoibiont  noi«*  conditions.  Sn , 
on  the  other  hand,  ie  superimposed  on  the  late  P  coda  and 
frequently  has  a  relatively  emergent  character  with  the  largest 
energy  on  the  vertical  component  often  associated  with  crustal 
multiples  (Kennett  1985).  Over  continental  areas  8n  usually 
propagates  well,  although  there  is  some  disruption  due  to  tectonic 
features  (liolnar  li  Oliver  1969). 

In  the  absence  of  widespread  digital  three-component  stations 
the  amplitude  comparison  of  Pn  and  Sn  is  normally  made  on  the 
vertical  component.  Thus  from  equations  (3.15),  the  vertical 
component  of  Pn  will  be  represented  by 
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with  a  phase  velocity  window  of  interest  from  7. 8-8. 3  km/s. 

For  shallow  sources  there  will  be  the  complication  that  the  time 
separation  of  direct  downward  radiation  and  surface  reflections 
will  be  too  short  to  allow  separate  identification  except  at  high 
frequency. 

From  (3.17,3.18)  the  vertical  component  of  the  Sn  phase  will 
be  described  by 
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for  phase  velocities  from  4.1— 4.7  km/s.  If  the  level  of 
heterogeneity  is  not  too  high  the  dominant  propagation 
remain  'in  plane  ,  with  a  consequent  simplification  of 
expression  for  the  vertical  component  of  Sn: 
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The  out-of-plane  scattering  terms  are  particularly  important  for 
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explosions,  «s  recognised  by  Gupta  &  Blandford  (1983),  since  in 
this  case  the  SH  energy  produced  at  the  source  will  be  low.  Even 
with  equal  efficiency  of  scattering  back  and  forth  between  SV  and 
SH,  there  will  be  transfer  of  energy  into  the  SH  field  until  the 
size  of  the  SV  and  SH  terms  tend  to  equalise.  For  earthquake 
sources  similar  cross-conversions  of  wavetype  will  occur  but  will 
be  harder  to  recognise  because  of  the  presence  of  significant 
SH  energy  generated  at  the  source. 

Mithin  the  mantle  the  character  of  the  propagation  terms  R,^  , 

PPP 

is  dictated  by  the  wavespeed  gradients  beneath  the  lioho. 

Even  a  slight  positive  gradient  is  sufficient  to  produce  an 
'interference*  head  wave  in  which  multiple  reflections  from  the 
gradient  constructively  interfere  to  give  a  larger  arrival  than 
is  possible  with  just  a  uniform  zone  (Henke  &  Richards  1980).  As 
the  range  increases  the  deepest  penetrating  'diving'  ray  separates 
from  the  bunch  of  multiple  reflections  near  the  crust-mantle 
interface  and  may  get  suppressed  by  interaction  with  a  wavespeed 
inversion  at  depth.  With  fine  scale  structure  in  'the  mantle,  the 
interaction  of  the  P  and  S  waves  will  be  different  because  of  the 
shorter  wavelengths  for  S. 

Comparing  (4.1)  and  (4.3),  we  see  that  they  have  a  similar 
structure  of  propagation  terms.  Loss  of  energy  by  out-of-plane 
scattering  will  tend  to  reduce  the  Sn  mantle  componment  somewhat 
and  also  the  transfer  to  the  vertical  component  is  less  efficient 
than  for  P  waves.  The  ratio  of  the  Pn  and  Sn  amplitudes  should 
however  be  dominated  by  the  ratio  of  the  source  terms  in  (4.1)  and 
(4.3). 

At  high  frequencies,  the  P  wave  contribution  to  (4.3) 
will  be  very  small  for  the  phase  velocities  in  Sn  (unless  the 
source  is  very  shallow),  so  that  the  Sn/Pn  ratio  becomes  a 
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comparison  of  S  and  P  radiation  at  similar  take  off  angles  from 
the  source.  At  lomer  frequencies,  and  for  shallow  sources,  the 
interference  of  PS  oonvereion  at  the  surface  with  the  direct 
downward  S  radiation  would  complicate  the  picture.  For 
explosions,  this  conversion  process  will  be  the  principal 
mechanism  for  the  genermtion  of  Sn. 

This  would  suggest  that  the  Pn/8n  ratio  would  perform  better 
as  a  discriminant  between  earthquakes  and  explosions  at  higher 
frequencies,  which  would  fit  in  with  a  number  of  observations 
(R.  Blandford  -  private  communication).  For  shallow  earthquakes 
the  discriminant  would  not  be  very  effective.  Clearly,  adequate 
azimuthal  coverage  is  needed  to  eliminate  complications  due  to 
nource  radiation  patterns. 

The  actual  behaviour  of  the  ratio  will  be  affected  by 
the  presence  of  heterogeneity.  Contamination  of  the  arrivals  by 
near — receiver  scattering  can  be  overcome  by  integrating  the 
envelope  of  the  waveforms  over  a  time  window  rather  than  just 
picking  a  maximum  value.  Where  available,  three  componcmt  data 
should  be  used  to  try  to  Stabilise  the  measure  of  the  ratio  of 
energy  in  Sn  and  Pn  against  distortion  due  to  energy  transfer 
between  the  vertical  and  horizontal  planes  due  to  scattering. 

Further,  in  order  to  equalise  the  sensitivity  of  the  Pn,Sn 
phases  to  structure  ( particularly  the  effects  of  heterogeneity)  it 
would  be  appropriate  to  work  at  a  common  wavelength  e.g.  make  a 
comparison  between  Pn  in  a  frequency  band  around  4  Hz  with  Sn  in  a 
frequency  band  around  2.5  Hz. 

4.2  Lg/Pn  4Mpiitud»  ratio 

Although  it  would  be  most  appropriate  to  compare  the  amplitudes 
of  the  two  crustal ly  guided  waves  Pq ,  Lg,  we  have  noted  above  that 
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the  Pg  phase  does  not  propagate  very  far  in  structures  with  high  P 
wave  velocity  at  the  surface  <such  as  shields).  For  utility  in  a 
wide  range  of  geologic  environments,  it  is  therefore  necessary  to 
use  the  Pn  phase  as  a  measure  of  the  strength  of  the  P  wave  field, 
except  at  short  ranges. 

With  a  comparable  approximation  to  (4.3)  the  vertical  component 
of  Lg  associated  with  in-plane  propagation  can  be  derived  from 
(3.12)  as 
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for  the  phase  velocity  range  3. 8-3.3  km/s.  The  propagation  term 
now  involves  reflection  from  the  crust/mantle  interface  in  R 

.4 

and  reverberations  within  the  waveguide  X .  The  coherence  of 
the  Lg  phase  depends  on  a  complex  pattern  of  constructive 
interference  which  is  normally  well  established  at  frequencies 
around  IHz,  except  where  some  major  barrier  blocks  propagation. 

Deep  reflection  soundings  of  the  crust  often  show  very  variable 
structure  at  the  crust/mantle  transition  along  a  profile  and  so  at 
high  frequency  (>  5  Hz)  the  propagation  effects  will  tend  to  get  more 
complex . 

Comparison  of  the  source  terms  in  (4.1)  and  (4.4)  shows  a 
similar  behaviour  to  that  for  Sn/Pn  with  a  simplification  at 
higher  frequencies,  although  now  the  Pn  waves  will  have  much  steeper 
take-off  angles  than  for  the  S  waves.  The  waves  comprising  Lg 
travel  at  angles  of  about  70  degrees  to  the  vertical,  whereas  the 
Pn  takeoff  angle  would  be  between  25  and  45  degrees  depending  on 
the  depth  of  the  source. 

The  Lg/Pn  amplitude  ratios  will  therefore  be  subject  to  two 
apposing  trends,  the  source  ratios  will  be  becoming  better  defined 
as  frequency  increases  but  the  propagation  complexity  will 
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increase  markedly  making  the  interpretation  of  the  ratio  rather 
difficult.  This  result  supports  the  discordant  observations  of  the 
value  of  the  amplitude  ratio  Lg/Pn  as  a  discriminant  noted  by 
Pomeroy  et  al  (1982). 

The  Lg  Mave  has  a  valuable  role  to  play  in  the  detection  of 
events  at  regional  ranges,  but  the  relatively  complicated 
propagation  characteristics  for  Lg  reduce  its  value  for 
discrimination . 

5.  Discussion 

In  this  paper  we  have  shown  how  an  economical  and  informative 
representation  of  the  seismic  wavefield  at  regional  distances  can 
be  made  in  terms  of  reflection  and  transmission  operators.  The 
expressions  we  have  derived  for  the  various  regional  phases  are 
valid  in  laterally  heterogeneous  and  anisotropic  media.  For  the 
particular  case  of  a  horizontally  stratified  media  they  can  be 
readily  evaluated  by  adapting  matrix  methods  in  the  transform 
domain  (see  e.g.  Kennett  (1983)).  In  the  presence  of 
heterogeneity  the  operators  can  be  simulated  by  allowing 
coupling  between  horizontal  wavenumbers  in  the  frequency  domain. 
Preliminary  calculations  show  that  it  is  possible  to  achieve 
significant  transverse  components  for  the  regional  S  wavefield 
from  an  explosive  source  as  has  been  frequently  observed 
(Blandford  1960),  but  which  is  not  predicted  by  a  stratified 
model . 

We  have  demonstrated  how  operator  representations  for 
particular  seismic  phases  can  be  applied  to  look  at  the 
theoretical  basis  for  discrimination  between  earthquakes  and 
underground  nuclear  explosions  using  regional  phases.  A  number 
of  empirical  discriminants  have  been  proposed  based  on  the 
characteristics  of  regional  phases  (see  Pomeroy  et  al  (1982)). 
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Our  phase  representations  give  a  way  of  examining  the  theoretical 
basis  of  such  discriminants  amd  may  allow  a  better  assessment  of 
their  limitations. 


Appendix:  Upper  and  Icjwer  crustal  heterogenei ty  zones 


In  section  2  we  have  shown  how  the  propagation  characteristics 
for  regional  phases  can  be  described  in  terms  of  the  regions 
above  and  below  the  source.  Frequently,  we  would  like  to  work 
with  fixed  regions  irrespective  of  source  depth,  so  that  we  can 
model  an  upper  crustal  zone  of  heterogeneity  down  to  a  level 

,  underlain  by  a  lower  crustal  zone  and  the  mantle.  Me  will 
assume  that  the  source  level  Zg  lies  above  and  will  show 
how  we  may  make  a  representation  of  the  wavefield  equivalent  to 
(2.11)  with  a  separation  of  the  propagation  characteristics  at 
the  level  k. 

If  we  had  a  source  at  the  level  K,  the  displacement  field  could 
be  written  in  an  analogous  form  to  (2.7) 
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C  I  -  R 
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and  we  could  then  make  a  separation  of  the  crustal  propagation 
terms  as  in  equations  (2.8-2.10)  in  terms  of  the  artificial  source 
level  K. 

In  order  to  incorporate  the  effects  of  the  true  source  at  S  we 
have  to  relate  the  apparent  source  terms  and  I/** 

to  the  actual  source  radiation.  Fortunately,  we  can  exploit  the 
formal  equivalence  with  reflection  matrix  results  for  splitting 
the  response  of  stratified  media  (Kennett  1983,  section  9.2).  Me 
find  for  the  net  downgoing  field 
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(A. 2) 


For  the  upgaing  Mevefield  et  K,  it  is  best  to  look  at  the  full 
propagation  term 
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3) 


which  we  may  recast  into  a  form  where  we  can  extract  the  direct 
upgoing  wave; 

H  to*  }  -  U  Tj*tV*  i 

+  fi  R*C  1  -  r!.*  R**  3*’  t/>‘  +  Ru  o'  ), 

”  •  “  *  -  u  ~ 

CIC 

In  general  we  would  not  expect  R^  to  be  too  large  so  that  the 

direct  wave  would  be  the  dominant  contribution.  In  both  (A. 2)  and 

(A. 4)  we  can  recognise  a  modified  source  term 

C  I  -  R^  R,  3  to*  ♦  R*  O*  > 

representing  the  effect  of  multiple  interactions  with 

heterogeneity  in  the  vicinity  of  the  source. 

On  the  receiver  side,  we  can  also  expand  out  the  upward 
^  He 

propagation  term  H  to  indicate  the  multiple  interactions 

with  the  shallow  heterogeneity 
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(A. 5) 


The  reverberation  term  will  be  responsible  for  most  of  the 
conversions  between  wavetypes  near  the  receiver. 

With  these  expressions  for  the  effect  of  shallow  heterogeneity, 
we  can  rewrite  the  operational  description  of  crustal ly  guided 
waves  (such  as  Pq ,  Lq )  into  the  form 

erujt  ~  ftc 
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with  the  main  propagation  characteristics  determined  by  the 

multiple  reverberations  within  the  crust  described  by 
*cc  4k  •i 

C  I  -  Rj  Ry  ]  •  For  those  phases  with  a  dominant  mantle  path 

(Pn,  Sn )  we  have  a  comparable  form 


s  K" 


KC  4K  kc 

C  I  —  R  R  1  T 


Hk  |r£ 


,  (A. 7) 

***  9 


with  the  mantle  propagation  operator  R^  modulated  by  possible 
crustal  reverberations,  and  the  possibility  of  strong  interactions 
with  heterogeneity  in  the  neai — source  (C**  )  and  neai — receiver 

J) 

(W  Ty*'  )  contributions . 
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Figure  Captions 

1.  Configuration  of  regions  and  surfaces  for  the  definition  of 
reflection  and  transmission  operators. 

2.  Schematic  representation  of  the  action  of  reflection  and 
transmission  operators  in  the  representation  of  the  seismic 
wavef ield. 

3.  Schematic  representation  of  the  operator  representation  for 
crustal ly  guided  waves  (3.2). 

4.  Schematic  representation  of  the  operator  representation  for 
the  main  mantle  phases  (3.14). 
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ABSTRACT 

The  Lg  wave  phase  which  is  of  considerable  interest  for  nuclear  discrimination 
problems  is  normally  observed  after  propagation  through  a  few  hundred  kilometres. 
This  phase  is  dominantly  guided  in  Ae  crustal  waveguide,  which  is  known  to  be  a 
region  with  very  significant  horizontal  variability  in  properties. 

The  effect  of  heterogeneous  crustal  structures  on  Lg  waves  has  been 
determined  by  using  a  "coupled-mode"  technique  in  which  the  local  seismic 
wavefield  in  the  real  medium  is  expressed  as  a  horizontally  varying  combination  of 
the  modal  eigenfunctions  of  a  stratified  reference  structure.  Departures  of  the  seismic 
properties  in  the  medium  from  those  of  the  reference  medium  lead  to  coupling 
between  the  various  amplitude  coefficients  in  the  modal  expansion.  The  evolution  of 
these  modal  weighting  factors  with  horizontal  position  are  described  by  a  coupled  set 
of  ordinary  differential  equations.  This  approach  provides  a  calculation  scheme  for 
studying  guided  wave  propagation  over  extended  distances,  at  frequencies  of  1  Hz 
and  above.  The  heterogeneity  models  which  have  been  us^  are  two-dimensional 
and  calculations  are  carried  out  fcK*  one  fiequency  at  a  dme. 

A  sequence  of  models  with  varying  levels  of  heterogeneity  have  been 
considered  in  order  to  determine  the  merits  and  limitations  of  the  computation 
scheme.  The  coupled  mode  technique  works  well  with  heterogeneous  models  in 
which  the  local  seismic  velocities  differ  from  the  stratified  reference  model  by  up  to 
2  per  cent  and  there  are  no  significant  distortions  of  the  main  discontinuities  (e.g.  the 
crust-mantle  boundary).  The  approach  can  be  used  for  higher  levels  of  heterogeneity 
and  with  distorted  interfaces  but  a  large  number  of  modes  needs  to  be  considered 
with  con^uent  high  computation  cc^ts.  If  the  level  of  heterogeneity  is  not  too  large 
then  the  interaction  between  modes  can  be  restricted,  rather  than  extending  over  £e 
whole  mode  set,  with  consequent  reduction  in  computation  cost. 

One  of  the  major  effects  of  crustal  heterogeneity  is  to  introduce  the  possibility 
of  smearing  out  the  main  amplitude  peak  in  the  Lg  wave  train  over  a  band  of  group 
velocities.  As  a  result,  an  effective  measure  of  the  energy  content  of  the  Lg  waves 
will  be  to  consider  the  integrated  amplitude  along  the  trace  beween  group  velocities 
of  3.6  and  3.3  km/s.  The  effects  of  heterogeneity  vary  between  different  parts  of  the 
Lg  wave  train  and  the  representation  of  the  wavefield  in  terms  of  modal 
contributions  allows  a  detail^  analysis  in  terms  of  the  group  velocity  components, 
which  can  be  illustrated  by  constructing  theoretical  seismograms  (with  a  narrow 
bandwidth  in  frequency)  for  the  heterogeneous  models. 
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INTRODUCTION 


Regional  S  phases  have  been  the  focus  of  considerable  attention  in  recent  years  in  the 
context  of  nuclear  discrimination  problems  (see  e.g  Pomeroy,  Best  &  McEvilly  (1982)) 
pc  Sn  and  Lg  wave  trains  are  well  established  by  about  200  km  from  the  source  and  can 
^ucntly  be  followed  out  to  ranges  of  1000  km  or  more.  These  S  phases  travel  through 
the  crust  and  uppermost  mantle  which  from  a  vari<^  of  studies,  are  known  to  be  regions 
with  very  considerable  horizontal  variability  in  properties.  In  order  to  improve 
understating  of  the  nature  of  regional  seismic  phases  and  the  way  in  which  the 
characteristics  imposed  by  the  source  may  be  modified  by  {xopagation  to  the  receiever,  we 
need  a  cmputadonal  procedure  which  will  allow  the  tracking  of  guided  seismic  wave 
propagatkm  through  a  horizontally  hetero^neous  crust  for  a  thousand  kilometres.  A 
suitable  candidate  is  the  coupled  mode  scheme  (tescribed  by  Kennett  (1984)  which  has 
already  been  used  with  some  success  in  understanding  the  characteristics  <rf  Le  wave 
behaviour  (Kennett  &  Mjdckeltveit  1984). 

For  regional  S  wave  trains,  a  representatkm  in  terms  of  a  limited  number  of  discrete 
modes  gives  a  economical  computatkmal  description  for  horizontally  stratified  media.  At 
each  frequency  only  a  limited  number  of  such  modes  need  be  considered.  In  a  laterally 
heterogeneous  medium,  at  a  single  frequency,  it  is  possible  to  represent  the  seismic 
displacement  And  traction  fields  widiin  the  varying  medium  as  a  sum  of  contributions  from 
the  modal  eigenfunctions  of  a  reference  structure,  with  coefficients  which  vary  with 
position.  The  evolution  of  diese  modal  expansion  coefficient  terms  with  horizontal  position 
can  be  described  by  a  set  of  coupled  partial  differential  equations  in  the  horizontal 
coordinates  (Kennett  1984).  The  cross-coupling  terms  between  different  modes  depend  on 
the  departures  of  the  heterogeneous  structure  from  the  stratified  reference  model.  These 
diffoences  are  not  required  to  be  very  small  but  must  not  be  such  as  to  completely  change 
the  nature  of  the  crustal  wave  guide.  Thus,  it  is  possible  to  accommodate  substantial 
localised  change  in  seismic  velocities  and  density,  but  more  difficult  to  allow  for  shifts  of 
more  Aan  2-3  km  in  the  position  of  major  interfaces  such  as  the  crust-mantle  boundary. 
Maupin  &  Kennett  (1987)  give  an  extei^ed  discussion  of  the  circumstances  in  which  the 
coupled  mode  approach  can  be  applied  to  two-dimensional  heterogeneous  models. 


STUDIES  OF  TWO-DIMENSIONAL  HETEROGENEITY  IN  THE  CRUST  AND 
MANTLE 

When  the  heterogeneity  within  the  medium  is  two-dimensional,  the  calculations  can  be 
recast  as  the  solution  of  non-linear  differential  equations  of  Ricatti  type  for  the  reflection 
and  ^smission  matrices  connecting  the  modal  expansion  coefficients  at  different 
positions  (Kennett  1984).  The  advantage  of  this  rearrangement  is  that,  for  each  frequency, 
the  boundary  conditions  on  the  differential  equations  are  simplified  for  a  generally 
heterogeneous  medium. 

If  we  adopt  a  reference  medium  which  docs  not  vary  with  horizontal  position,  the 
individual  mode  contributions  propagate  independently  in  that  structure.  However,  once 
the  imrperties  of  the  true  medium  differ  from  the  reference,  the  independence  of  modal 
propagation  is  lost  The  effect  of  the  heterogeneity  enters  into  the  differendal  equations  for 
the  modal  expansion  coefficients  via  a  coupling  matrix  whose  dimensions  are  dictated  by 
the  number  of  modes  included  in  the  calculation.  The  choice  of  the  number  of  modes  is  of 
coiuiderable  importance.  All  the  significant  wavenumber  conqxMients  for  the  seismic  phase 
of  interest  should  be  included,  as  well  as  an  allowance  for  steeper  angles  of  propagation 
than  would  be  present  in  the  reference  medium,  in  order  to  allow  for  scattering  effects.  The 
demands  of  computation  are  that  the  number  of  modes  should  be  kept  as  small  as  possible, 
since  the  computation  time  depends  on  the  square  of  the  number  of  modes.  On  the  other 
hand  if  the  truncation  is  too  tight,  a  poor  representation  of  the  displacement  field  can  ensue 
and  the  results  are  of  limited  use. 

In  figure  1  we  illustrate  the  reference  model  (ARANDA)  used  for  most  of  the  Lg  mode 
calculations,  and  the  modal  eigenfunctions  for  the  22  Love  and  Rayleigh  modes  included  in 
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calculations  at  1.5  Hz.  'Hie  vertical  component  of  the  Rayleigh  modes  are  displayed  to 
illustrate  their  very  similar  depth  dependence  to  the  horizont^  component  of  the  Love 
modes.  The  cut-off  criterion  was  based  on  the  effective  depth  penetration  of  the  modes. 
The  fundamental  and  first  higher  modes  are  confined  to  the  sedimentary  zone.  Modes  2  to 
IS  are  the  main  contributors  to  the  Lg  phase  with  nearly  all  their  energy  trapped  in  the 
crust.  Mode  16  begins  to  have  significant  displacement  in  the  mantle,  and  the  remaining 
modes  carry  their  main  energy  in  ^e  mantle  and  represent  the  Sn  phase.  The  truncation  is 
made  at  a  phase  velocity  of  4.54  km/s,  corresponding  in  this  model  to  a  penetration  depth 
of  about  100  km.  Features  in  the  wave  propagation  process  corresponding  to  higher  phase 
velocities  than  4.54  knVs  cannot  be  represented  by  these  coupled  noiode  calculations. 

In  order  to  give  a  good  account  of  the  guid^  wave  propagation,  we  must  make  sure 
that  all  the  signiHcant  wavenumber  components  for  the  seismic  phases  of  interest  are 
included,  as  well  as  making  an  additional  allowance  for  steeper  angles  of  propagation  in 
order  to  allow  for  scattering  effects.  This  means  that  we  cannot  just  use  those  modes 
which  have  the  bulk  of  their  energy  within  the  crust  but  must  also  include  waves  with 
higher  phase  velocity  whose  energy  in  the  stratified  reference  model  lies  dominantly  in  the 
upper  mantle.  At  1.5  Hz,  we  need  to  carry  at  least  the  22  modes  illustrated  if  we  are  to 
allow  for  wave  interaction  processes  generating  waves  travelling  at  angles  up  to  40*  from 
the  horizontal  in  the  midcrust  and  even  then  we  will  not  give  a  full  account  of  waves 
travelling  at  steeper  angles  to  the  horizontal. 

Significant  perturbations  of  the  major  interfaces  in  the  model  lead  to  extensive  coupling 
between  modes  and  can  induce  reladvely  steep  angles  of  propagation  (Maupin  &  Kennett 
1987).  If  such  features  are  required  in  the  nnodels  to  be  conside^,  then  a  broad  sweep  of 
modes  must  be  taken  for  the  reference  model  (for  the  Lg  case  at  1.5  Hz,  at  least  30).  With 
the  benefit  of  hindsight,  it  would  appear  that  the  structural  model  used  in  the  calculations  of 
Kennett  &  Mykkeltveit  (1984)  on  the  influence  of  graben  structures  on  Lg,  was  as  large  a 
deviation  fiom  the  reference  rnodel  as  could  be  treated  effectively  using  their  20  mode  set. 

The  coupled  mode  representation  is  well  suited  to  heterogeneous  models  where  the 
behaviour  consists  of  relatively  random  variations  about  the  properties  of  the  reference 
model  (see  e.g.  fig  2).  In  this  case,  a  single  incident  nuxle  will  interact  with  a  limited 
number  of  neighbouring  modes  so  that  the  coupling  matrix  is  dominated  by  the  diagonal 
elements,  with  a  limited  effective  bandwidth.  With  the  aid  of  the  coupled  mode  technique, 
a  study  has  been  made  of  the  way  in  which  the  modal  field  for  Lg  waves  is  affected  by 
various  levels  of  distributed  heterogeneity  for  frequencies  up  to  2Hz  over  paths  of  1(X)  to 
10()0  km  long  with  both  long  and  short  horizontal  scales  of  heterogeneity  within  the  crust 
and  mantle. 

Broad  Scale  Heterogeneity 

We  will  firstly  describe  simulations  of  long  range  propagation,  for  which  a  sequence 
of  different  heterogeneity  models  were  constructed  by  sp^ifying  the  seismic  velocities  at 
40  km  horizontal  intervals,  for  a  sequence  of  depths  in  a  vertical  section  and  then  using 
bicubic  spline  interpolation  within  each  layer.  For  the  distributed  heterogeneity  models  the 
velocity  values  were  generated  with  a  random  perturbation  to  the  reference  v^ue  within  a 
prescribed  range  of  variation.  The  velocity  vdues  were  then  smoothed  horizontally  by 
applying  a  moving  average  over  3  points  to  avoid  extreme  variations  with  short  horizontal 
scales.  The  horizontal  smoothing  process  imposes  a  typical  horizontal  scale  length  of 
around  100  km,  and  the  vertical  scale  varies  wi^in  the  niodel  becoming  l^er  in  the  upper 
mantle.  The  velocity  distribution  is  arranged  so  that  there  are  no  deviations  from  the 
reference  structure  at  the  ends  of  the  model  (0  and  1(X)0  km).  This  enables  a  direct  physical 
interpretation  of  the  mode  coupling  results,  after  passage  through  the  heterogeneous  region, 
in  terms  of  the  reflection  and  transmission  of  the  modes  of  the  reference  structure. 

In  figure  2  we  illustrate  a  set  of  models  (D,  E,  F)  with  increasing  crustal  heterogeneity. 
Model  D  was  constructed  by  imposing  a  perturbation  to  the  velocities  at  the  specific  knots 
of  the  interpolation,  chosen  from  a  uniform  random  distribution  on  ±  2  p^ent.  After 
horizontal  smoothing  by  averaging  over  3  horizontal  knots,  the  resulting  perturbations  in 
model  D  have  about  ±1  per  cent  heterogeneity  in  the  crust  with  a  horizontal  scale  length  of 
about  1(X)  km,  and  ±1  per  cent  heterogeneity  has  been  introduced  in  the  mantle  as  well. 
Models  E  and  F  retain  the  same  mantle  heterogeneity  structure  as  model  D  but  have  higher 


cnisul  heterogeneity  amplitudes,  in  the  initial  perturbation  stage  up  to  ±5  per  cent  variation 
from  the  reference  model  was  allowed  for  model  E  and  ±10  per  cent  for  nxxkl  F.  After 
horizontal  smoothing,  die  variations  from  the  reference  model  are  of  the  order  of  ±2  per 
cent  for  model  E  and  ±5  per  cent  fm'  model  F.  We  shall  also  consider  noodel  C  with  the 
same  crustal  structure  as  D,  but  no  heterogeneity  below  the  crust-mantle  boundary. 

After  propagating  a  large  horizontal  distance  in  a  heterogeneous  model  ^e  energy 
originally  in  a  single  incidoit  mode  will  no  longer  be  confined  to  that  mode,  and  indeed 
some  energy  iruiy  be  r^ected  back  by  the  heterogeneity,  with  a^ain  the  possibility  of 
conversion  between  modes.  When  the  level  (tf  heterogeneity  is  relatively  low  (of  the  order 
(tf  1  per  cent  deviations  from  the  stratified  reference  model  -  as  in  model  D)  transmission 
effects  dominate  and  the  level  of  reflected  energy  is  negligible.  Even  with  rather  larger 
levels  (tfbeterc^eneity  die  reflection  from  the  wh^  model  is  small  but,  as  ptmted  out  by 
Maupin  &  Kennett  (1987),  it  is  necessary  to  retain  apparmdy  reflected  waves  in  the  course 
of  the  computadon  in  or^  to  get  an  accurate  calculadon  of  transmission  effects.  The 
Rkatd  equations  are  much  sinqiler  when  reflecdon  can  be  neglected  and  so  computadon 
dme  can  be  reduced. 

PCX'  any  pardcular  heterogeneous  structure,  it  is  therefore  necessary  to  undertake  a 
preliminary  t^culadon  including  both  reflecdon  and  transmission  effects  befexe  it  can  be 
assumed  dut  reflecdon  can  be  ignored.  The  control  of  numerical  accuracy  in  the  course  of 
the  calculations  is  rather  difficult,  although  for  perfeedy  elastic  models  there  is  the 
possibility  of  monitexing  the  constancy  of  the  total  energy  in  reflection  and  transmission 
associated  with  an  individual  incident  nxxie.  To  preserve  full  numerical  precision  quite 
small  steps  have  to  be  taken  in  the  hcxizcxital  dir^on:  at  1.5  Hz  the  step  should  be  no 
larger  dum  0.5  bn  if  bodi  reflected  and  transmitted  waves  are  considered.  For  a  number  of 
mc^ls  a  somewhat  larger  step  length  gave  reasonable  results  for  the  case  of  transmission 
akme. 

The  numerical  calculations  presented  in  this  paper  have  been  performed  with  the 
inclusion  of  both  reflecdon  and  transmission  effects,  and  energy  conservation  was  satisfied 
to  better  than  0.25  per  cent  over  KXK)  bn  of  propagation  through  heterogeneous  structure. 
To  reduce  possible  ernxs  due  to  truncation  of  the  modal  sum,  the  calculations  were  carried 
out  using  a  modal  suite  at  least  5  modes  larger  than  requir^  by  the  truncation  criterion. 
The  reflection  and  transmission  matrices  for  the  full  heterogeneous  region  were  then 
trancated  before  di^lay  or  odier  analysis. 

For  an  incident  mode  at  0  bn  the  effect  of  KXX)  km  of  heterogeneity  is  that  the  energy 
originally  in  the  mode  is  no  longer  confined  to  that  particular  mode.  In  figure  3  we  show 
the  transmission  matrix  for  both  Love  and  Rayleigh  modes  with  the  22  mode  set 
appropriate  to  the  ARANDA  reference  model  at  1.5  Hz,  after  passage  through  the 
heterogeneity  model  D.  Each  column  of  the  matrix  corresponds  to  the  modal  coefficients 
produced  by  the  incidence  of  a  single  mode  of  amplitude  100  at  0  km.  We  see  that  the 
behaviour  is  dominated  by  the  diagonal  elements,  but  that  significant  mode  conversions 
occur  producing  sizeable  off-diagonal  elements.  For  the  Lg  type  modes  (2-15)  the 
bandwidth  for  signifleant  interaction  is  typically  ±  2  modes.  The  strong  interaction  of 
modes  2  and  3  arises  because  of  their  similar  shape  in  the  top  10  km.  For  the  Sn  modes 
(16-21),  the  eigenfunctions  are  very  similar  and  there  is  strong  interaction  due  to  mantle 
perturbation  extending  to  200  km  (the  lower  portion  of  the  heterogeneitv  models  are  not 
shown  in  flg  2).  However  the  separation  fnxn  the  Lg  modes  is  suibng;  thore  is  very  little 
interaction  except  fex*  modes  15  and  16  which  share  some  of  the  characteristics  of  each 
group.  Tlie  behaviour  for  Love  and  Rayleigh  waves  are  similar,  although  they  display 
some  difference  in  their  sensitivity  to  heterogeneity  because  the  coupling  terms  in  die  Ricatti 
equations  depend  on  different  combinations  of  physical  parameters. 

As  the  fRX]uency  increases  the  number  of  modes  ne^  to  cover  the  same  phase  velocity 
range  increases,  which  causes  some  problems  in  combining  the  results  from  different 
frequencies.  The  increased  size  of  the  differential  equation  system  also  means  an  increase 
in  computation  time  as  the  square  of  the  number  of  modes  if  all  mode  interactions  have  to 
be  considered.  Ftxtunately,  for  the  same  model,  the  bandwidth  of  interaction  increases 
only  slowly  with  frequency.  For  Lg  waves  propagating  through  models  similar  to  D  with 
about  1  per  cent  heterogeneity,  computations  can  be  restricted  to  a  band  of  7  modes  either 
side  of  the  target  mode  without  appreciable  error  for  frequencies  up  to  2  Hz.  At  1 .5  Hz, 
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with  22  modes  included,  this  bandwidth  restriction  has  the  effect  of  reducing  computation 
time  by  nearly  a  factor  of  three. 

The  required  bandwidth  increases  with  the  level  of  heterogeneity  in  the  iiKxlel  (see 
figure  5).  For  model  D  in  figure  2,  a  bandwidth  of  ±6  modes  around  the  target  mode  is 
sufficient  for  both  Lg  and  Sn  modes.  However  to  account  for  the  interaction  between 
modes  for  model  E  it  is  desirable  to  have  at  least  a  bandwidth  of  ±8  modes  included  in  the 
calculation.  For  the  highest  level  of  heterogeneity  illustrated  in  fig  2  (model  F)  a  bandwidth 
of  ±10  modes  is  definitely  insufficient  for  full  accuracy.  For  very  heterogeneous  models 
there  is  also  the  possibility  of  high  angle  propagation  effects  and  so  the  size  of  die  mode  set 
has  to  be  increa^.  Thus  although  the  coupled  mode  approach  is  not  confined  to  low  level 
heterogeneity,  the  method  is  most  effecdve  when  the  structure  does  not  deviate  too  far  from 
the  reference.  A  convenient  woridng  limit  would  seem  to  be  about  2  per  cent  deviation  as 
in  model  E. 

In  order  to  understand  the  effect  of  propagation  through  a  heterogeneous  model  on  the 
nature  of  regional  phases  we  need  to  understand  the  characteristics  of  the  surface  wave 
modes  in  the  reference  model.  In  figure  4  we  display  the  group  slowness  behaviour  as  a 
function  of  frequency  up  to  2  Hz  for  the  first  40  Rayleigh  modes  on  the  ARANDA 
reference  structure.  Ilie  advantage  of  displaying  the  group  slowness  (the  reciprocal  of  the 
group  velocity)  is  that  the  various  arrivals  can  be  recognised  in  the  time  relationship  that 
Aey  have  on  a  seismic  record 

The  Sn  phase  can  be  recognised  in  figure  4  as  the  superposition  of  many  slowness 
maxima  and  minima  at  a  slowness  around  0.22  sflan  (4.S  kin/s).  All  40  modes  contribute 
to  this  tightly  defmed  band  of  slownesses  which  clearly  separates  from  the  other  classes  of 
arrival.  The  onset  of  the  Lg  phase  at  a  slowness  of  0.286  s/km  (3.S  km/s)  is  associated 
with  the  superposition  of  a  number  of  relatively  broad  slowness  minima  arising  from 
relatively  low  order  modes.  Following  this  onset  is  a  more  tangled  skein  of  maxima  and 
minima  with  slownesses  less  than  0.3(M  s/km  (i.e.  group  velocities  greater  than  3.3  km/s), 
associated  with  about  four  modes  at  each  fr^uency,  which  will  provide  the  main  Lg 
arrival.  Somewhat  later,  we  have  a  sequence  of  well  defined  slowness  maxima  associated 
with  each  mode  in  turn  as  the  frequency  increases.  These  maxima  represent  the  Airy  phase 
for  each  mode  and  will  be  significant  contributors  to  the  Lg  coda;  they  arise  physically  from 
multiple  reflections  at  angles  to  the  vertical  just  greater  than  the  critical  angle  for  the 
crust-mantle  boundary. 

From  the  modal  eigenfunction  patterns  displayed  in  figure  1,  we  can  recognise  that 
those  modes  (3-6)  which  contribute  most  to  the  onset  of  Lg  at  1.5  Hz  have  the  bulk  of  their 
energy  confined  to  the  middle  and  upper  crust.  The  main  Lg  arrivals  come  from  iiKxles 
7-10  which  sample  the  whole  crust.  The  Airy  phases  for  modes  14,  15  arise  just  before 
the  transition  to  energy  transport  in  the  mantle. 

The  effect  of  varying  the  heterogeneity  level  in  the  structural  models  can  be  well 
illustrated  by  comparing  the  energy  ^stribution  across  the  modal  sequence  for  a  single 
incident  mode.  In  figure  5  we  show  this  energy  distribution  for  a  set  of  incident  nnodes  at  0 
km,  at  a  frequency  of  1.5  Hz,  chosen  to  represent  different  parts  of  the  Lg  wave  phase. 
We  consider  horizontal  transmission  through  1000  km  of  structure  for  four  different 
models  C,  D,  E  and  F  with  increasing  levels  of  heterogeneity.  Model  C  has  1  per  cent  of 
hetrogeneity  confined  to  the  crust.  M^els  D,  E  and  F  are  illustrated  in  figure  2,  and  have 
the  same  heterogeneity  model  in  the  mantle  with  about  1  per  cent  variation  from  the 
reference  model  (ARANDA,  fig  1),  but  increasing  levels  of  crustal  heterogeneity  (±1  per 
cent  for  D,  ±2  per  cent  for  E,  and  ±5  per  cent  for  F). 

At  1.5  Hz,  mode  4  is  a  principal  contributor  to  the  onset  of  the  Lg  wave  train.  From  fig 
5a  we  see  that  for  the  incidence  of  this  single  mode  at  0  km,  the  majority  of  the  energy  is 
carried  in  the  original  mode  over  the  full  1000  km  propagation  path  for  the  structures  with 
up  to  2  per  cent  variation  away  from  the  reference  model  (C,D,E).  However,  the 
proportion  of  energy  in  the  original  mode  decreases  as  the  crustal  heterogeneity  increases. 
As  expected  the  mantle  heterogeneity  has  very  little  influence  on  this  group  of  interacting 
modes  which  is  principally  confin^  to  the  crust.  Once  we  reach  the  highest  level  of 
heterogeneity  (n^el  F)  the  pattern  of  the  energy  distribution  is  markedly  changed,  the 
energy  maximum  is  shifted  over  two  modes  to  mo^  6  which  has  a  slightly  different  group 
velocity  and  energy  is  spread  widely  across  the  whole  suite  of  crustal  modes.  There  is  very 
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Ikde  intencdon  with  the  fundamental  or  first  higher  inodes  because,  as  can  be  seen  from 
figve  U  their  eno-gy  is  confined  to  the  near  surface  and  so  is  only  sensitive  to  a  small  part 
of  die  total  cnistal  heieiogendty. 

For  a  mode  in  the  main  Lg  wave  arrivals  (mode  9),  we  see  from  fig  5b  that  for  the  two 
modds  widi  relatively  modest  levels  of  hetoogeneity  (CD)  the  energy  is  cmcentrated  over 
teee  aeighboming  modes  centered  on  the  incident  mode.  These  dim  modes  have  group 
vdodties  varying  by  about  O.IS  knVs  and  so  the  effect  of  die  this  level  of  heterogeneity 
wiD  be  to  give  a  more  ddiuse  maximum  in  the  wave  amplitude.  The  ]»esence  mande 
heienigendty  now  Insadight  effect  since  it  ispossiUe  to  couple  into  modes  with  some 
energy  penetration  into  the  uppermost  mantle.  For  die  higher  levek  of  crustal  betpogeneity 
the  incKlent  energy  gets  out  over  a  broader  range  of  modes  (6  are  significant  for 
model  E  and  at  feast  10  ^  model  F)  and  the  effect  on  the  Lg  wavetrain  will  be  more 
profound. 

As  a  result  of  the  di^iersal  energy  across  a  number  (tf  modes,  for  modest  levels  of 
hetnogeneity  (around  1  per  cent),  the  best  measure  of  Lg  wave  strength  will  be  the 
integn^  energy  in  a  time  window  ^panning  group  velocities  of  about  3.6-3.3  km/s  for  a 
typ^  continental  situation;  since  die  modes  with  diese  properties  have  their  dominant 
interactions  within  the  group. 

Similar  results  arise  for  the  Lg  wave  coda  (incident  mode  14)  though  this  case,  where 
the  incident  mode  is  most  sensitive  to  lower  crustal  structure,  is  not  as  stron^y  affected  by 
the  heterogeneity  as  for  mode  9.  Even  so,  as  we  see  in  figure  Sc,  with  increasing 
heterogeneity  levels  the  representation  of  the  propagation  via  the  modes  of  the  reference 
structure  requires  extoisive  mode  coupling.  Fch*  the  higher  heterogeneity  levels  there  is 
significant  energy  shift  into  Sn  type  modes,  which  we  can  explain  physictdly  as  occurring 
from  local  Ganges  in  the  critical  angle  at  the  crust-mande  boundary.  Modes  like  14  travel 
close  to  the  critical  angle  in  the  reference  model  and  so  a  slight  change  of  conditions  can 
lead  easily  to  energy  transmission  into  the  mantle,  with  a  consequent  shift  of  amplitude  to 
modes  like  17  arid  18  (in  the  Sn  suite).  Energy  is  conserved  across  the  whole  modal  field 
but  is  ledistribuied  between  the  contributions  of  individual  modes. 

Fn-  incident  modes  4  and  9,  the  nature  of  the  eno-gy  spread  induced  across  the  modes 
of  the  reference  stracture  has  bom  principally  to  couple  Lg  type  modes  together  and  so  the 
pattern  of  energy  spread  shows  a  skew  towards  higher  mode  numbers  for  mode  4  and 
lower  numbers  fOT  mode  9.  Whereas,  for  incident  mode  14,  the  coupling  is  largely  into  a 
the  imnvidiaftR  neighbouring  modes  with  similar  character  and  also  into  the  Sn  nv^es. 

We  should  note  that  even  for  the  moderately  heterogeneous  model  D  (fig  3)  there  has 
been  a  significant  modification  of  the  distribution  of  energy  between  nxxles  after  passage 
tluxHigh  1000  km  of  structure.  Such  modifications  over  a  band  of  frequencies  will  lead  to 
theoretical  seismograms  which  will  differ  significantly  from  the  predictions  for  the 
hoiixontally  stratified  reference  model. 

We  can  simulate  the  effect  of  heterogeneity  on  theoretical  seismograms,  by  modulating 
the  modal  coefficients  by  the  appropriate  transmission  matrix.  In  order  to  reduce 
compuutional  effort,  whilst  displaying  a  useful  result,  we  have  considered  seismograms 
for  a  narrow  frequency  band  and  interpolated  the  modal  transmission  matrix  between 
calculations  for  just  a  few  frequencies.  We  restrict  attention  to  explosive  sources  to  simplify 
the  source  behaviour.  In  figure  6,  we  show  seismograms  at  1050  km  for  a  source  at  5km 
depth  in  the  ARA^A  model  with  a  centre  frequency  of  1  Hz.  The  upper  trace  is  for  the 
htxizontally  stratified  reference  model  and  is  dominated  by  Lg  waves  with  relatively  weak 
Sn.  The  two  lower  traces  include  the  effect  of  transmission  through  10(X)  km  of  the 
heterogeneous  models  D  and  E.  The  seismograms  are  normalised  to  the  same  peak 
amplioi^  and  relative  amplitudes  are  indicated  to  the  right  of  figure  6.  The  behaviour  of 
these  theoretical  seismograms  illustrates  clearly  the  character  of  the  intermode  interactions 
we  have  already  ^scussed.  The  ARANDA  model  has  a  weak  gradient  for  S  velocity  in  the 
mantle  and  so  the  Sn  phase  is  relatively  weak  fOT  the  horizontally  stratified  nKxlel. 
However  the  coupling  between  the  Lg  and  Sn  wavegudes  induced  by  the  heterogeneity 
leads  to  an  increase  in  Sn  amplitude,  principally  at  the  expense  of  the  Lg  coda.  Whilst 
model  D  with  only  1  per  cent  hetrogeneity  shows  most  change  in  the  onset  of  Lg,  the 
higher  level  of  heterogeneity  in  model  E  has  lead  to  a  smearing  out  of  the  most  energetic 
part  of  the  Lg  phase  with  an  amplitude  reduction  of  about  10  per  cent.  The  Lg  wavetrain 
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can  be  viewed  either  as  a  superposition  of  many  multiple  reflections  within  the  crustal 
waveguide  or  as  a  interference  of  many  modes.  The  effect  of  heterogeneity  is  to  change  the 
relative  arrival  times  of  the  multiple  reflections  and  thus  the  phases  of  the  different  modes 
so  that  the  interfence  pattern  changes.  Because  there  is  a  delicate  balance  between 
constructive  and  destructive  interference,  relatively  small  changes  in  timing  can  have  a 
substantial  effect  on  the  waveform. 

These  results  show  that  propagation  over  hundreds  of  kilometres  through  structures 
with  relatively  long  wavelength  of  heterogeneity  can  have  substantial  effect  on  the  character 
of  the  Lg  wave  train.  However,  in  addition  to  large  scale  effects,  there  is  abundant 
evidence  for  small  scale  heterogeneity  on  a  scale  of  tens  of  kilometres  with  amplitudes  of  a 
few  percent  (see  e.g.  Aki  1981).  We  need  to  know  what  is  the  effect  of  such  small  scale 
heterogeneity  and  how  it  may  affect  long  range  propagation. 

Small  scale  heterogeneity 

In  order  to  look  at  concentrated  local  heterogeneity  we  have  considered  the  effects  of 
the  same  velocity  perturbations  as  before  but  on  a  compressed  distance  scale.  We  have 
taken  the  same  velocity  values  as  in  models  D  and  E  above  (fig  2)  but  reduced  the 
horizontal  scale  by  a  factor  of  10  so  that  in  the  new  models  H  (approximately  ±1  per  cent 
variation)  and  I  (approximately  ±  2  per  cent)  the  heterogeneity  is  concentrated  in  a  100  km 
zone.  The  velocity  sampling  is  now  at  4km  horizontal  intervals.  This  yields  velocity 
perturbations  of  a  few  percent  with  a  horizontal  heterogeneity  scale  of  around  10  km, 
representing  the  smaller  scale  features  commonly  found  in  the  crust. 

For  these  scales  of  heterogeneity,  as  frequency  increases  and  consequently  wavelength 
diminishes,  the  size  of  the  off-diagonal  elements  in  the  modal  coupling  matrix  increase  at 
the  expense  of  the  diagonal  elements  which  will  induce  a  broader  spread  in  the  group 
slownesses  associated  with  the  propagation. 

The  heterogeneity  model  H  has  relatively  little  effect  on  the  modes;  there  is  a  spread  of 
ener^  away  from  the  incident  mode  after  a  100  km  passage  through  the  nnodel,  but  this  is 
principally  to  nearest  neighbour  modes  which  will,  in  general,  have  similar  character.  Once 
the  heterogeneity  level  increases  (to  around  ±  2  per  cent  in  model  I)  the  behaviour  is  not  so 
simple .  The  transmission  matrix  is  sdll  diagcxial  dominated  but  there  is  now  coupling  over 
a  span  of  a  number  of  modes.  These  effects  can  be  illustrated  well  by  once  again  forming 
theoretical  seismograms  with  a  narrow  frequency  band.  To  allow  direct  comparison  with 
figure  6  we  have  once  again  considered  the  seismograms  at  a  distance  of  10^  km  with  a 
centre  frequency  of  1  Hz. 

In  figure  7  we  show  the  theoretical  seismograms  at  1050  km  including  the  effect  of 
passage  through  1(X)  km  of  small  scale  heterogeneity;  the  upper  trace  is  far  the  hraizontally 
stratified  reference  model.  For  model  H  (±  1  per  cent)  the  difference  in  waveform  is  slight 
and  can  only  be  seen  on  careful  inspection  but  the  amplitude  has  decreased  by  nearly  5  per 
cent.  However,  with  ±  2  per  cent  heterogeneity  as  in  nnodel  I  there  is  a  noticeable  change 
in  waveform  and  a  small,  but  significant,  transfer  of  energy  ahead  of  the  main  Lg  phase. 
The  amplitude  maximum  is  now  more  sharply  defined  and  large  than  in  the  absence  of 
heterogeneity  due  to  favourable  constructive  interference.  The  coupling  between  modes  is 
at  its  strongest  for  the  lower  orders  and  this  is  reflected  by  the  change  in  the  nature  of  the 
onset  of  Lg.  Thus,  a  substantial  change  in  the  character  of  the  Lg  wavetrain  can  be 
produced  by  only  1(X)  km  of  moderate,  small-scale,  heterogeneity.  So  that  it  is  clear  that 
the  cumulative  effect  of  multiple  scattering  within  the  crust  can  be  substantial  and  that 
observed  waveforms  can  be  heavily  shaped  by  the  propagation  path  as  well  as  the  source 
characteristics  (even  in  the  absence  of  major  lateral  heterogeneity). 

DISCUSSION 

The  results  of  the  studies  of  Lg  wave  propagation  in  heterogeneous  models  for  the 
crust  and  upper  mantle  show  that  reflection  from  distributed  heterogeneity  is  not  very 
important.  However,  both  reflection  and  transmission  processes  should  be  included,  if  at 
all  possible,  in  the  numerical  calculations  in  order  to  ensure  the  accuracy  of  estimates  of 
modal  transmission.  The  coupled  mode  technique  works  well  with  heterogeneous  models 
which  differ  from  a  stratified  reference  model  by  up  to  2  per  cent  with  a  limited  bandwidth 
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oi  interaction  between  inodes.  The  i^roach  can  be  used  for  higher  levels  of  heterogeneity 
but  a  large  numbor  of  inodes  need  to  be  considered  with  consequent  high  computation 
costs. 

One  of  the  major  effects  of  crustal  heterogeneity  is  to  introduce  the  possibility  of 
smeumg  out  the  main  amfriiuide  peidc  in  the  Lg  wave  train  over  a  band  of  group  velocities. 
As  a  rewk.  an  ^ecttve  measae  of  the  energy  conteni  of  the  Lg  waves  will  be  to  consider 
die  imegnted  ami^tude  akmg  die  ttaces  beween  group  velocities  of  3.6  and  l3  kWs.  The 
influeaoe  of  maiMle  hyityncity  on  Lg  wave  prop^atimi  was  found  to  be  quite  small. 
However,  any  convcraon  inm  Sn  type  modes  will  give  arrivals  which  will  have  iqipan^nt 
group  veloddes  of  between  4.3  and  kn^,  d^iending  on  the  position  where  conversion 
occurs.  Sndi  arrivals  will  appear  in  a  portion  of  die  seismic  record  which  would  be 
jxedicted  to  be  very  qmet  in  stratified  medium  calculations,  and  as  a  result  may  appear  to  be 
more  prondneot  than  expected.  The  increase  of  energy  in  this  time  interval  for  the 
h^eiogeneous  models  goa  some  way  towards  leamciling  theoretical  seismograms  with 
the  diaracto’of  obsov^  regional  phaxs. 
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HGURE  CAPTIONS 

Figure  1:  ARANDA  reference  model  with  Love  and  Rayleigh  mode  eigoifuncticHis  for 
phase  velocity  less  than  4.54  km/s,  at  frequency  1.5 

Figtro  2:  Contour  plots  of  the  s^uence  of  heterogeneous  velocity  models  used  in  the 
studies  of  Lg  wave  propagation.  The  heavier  lines  indicate  positive  perturbations  away 
fiom  the  AI^NDA  reference  model,  and  the  contour  intoval  is  0.5  per  cent 
perturbation. 

Figure  3:  Reptesentadon  of  the  transmissicm  matrix  for  Love  and  Rayleigh  modes  at  1.5  Hz 
after  pa^ge  duough  1000  km  of  model  D.  Anqilitude  values  less  tl^  1  are  left  blank. 

Figure  4:  Group  slowness  beviour  as  a  function  of  frequency  for  the  first  40  Rayleigh 
noodes  of  the  ARANDA  reference  model  illustrated  in  figure  1. 

Figure  5:  Energy  distribution  across  die  Love  and  Rayleigh  modes  for  the  ARANDA 
r^erence  model  at  1.5  Hz,  after  transmissitHi  thnwgh  1000  km  of  heterogeneous 
structure.  The  symbols  denote  die  different  hetrogeneity  models,  as  discuss^  in  the 
text:  a)  incident  mode  4,  b)  incident  mode  9,  c)  incident  noode  14. 

Figure  6:  Theoretical  seismograms  at  1050  Ian  for  a  narrow  ftequency  band  around  1  Hz 
illustrating  the  effect  of  teoad  scale  heterogeneity  on  die  regimial  S  wavetrain.  The 
uiqier  trace  is  for  die  horizontally  stratified  ARANDA  reference  model  and  the  two 
lower  traces  fOT  propagation  thrwgh  1000  km  of  heterogeneiQr  models  D  (±  1  per  cent) 
and  E  (±  2  per  cent).  The  seismc^iams  are  normalised  to  the  same  peak  an^litude  and 
the  relative  an^litudes  are  indicated  to  die  ri^t  of  each  seismogram. 

Figure  7:  Thecaedcal  seismograms  at  1050  km  for  a  narrow  frequency  band  around  1  Hz 
illustrating  the  effect  sr^  scale  heterogeneity  on  the  regicmal  S  wavetrain.  The 

upper  trace  is  fix'  the  hcdzuiitally  stratified  AR/d4DA  reference  naodel  and  die  two 
lower  traces  for  propagadon  through  100km  of  heterogeneity  models  H  (±  1  per  cent) 
and  I  (±  2  per  cent).  The  seismograms  are  normalised  to  the  same  peak  ampli^e  and 
the  reladve  amplitudes  are  indicated  to  die  right  of  each  seismogram. 
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